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Abstract 

We describe a large class of chemical reaction networks, those en- 
dowed with a subtle structural property called concordance. We show 
that the class of concordant networks coincides precisely with the class 
of networks which, when taken with any weakly monotonic kinetics, 
invariably give rise to kinetic systems that are injective — a quality 
that, among other things, precludes the possibility of switch-like tran- 
sitions between distinct positive steady states. We also provide persis- 
tence characteristics of concordant networks, instability implications 
of discordance, and consequences of stronger variants of concordance. 
Some of our results are in the spirit of recent ones by Banaji and 
Craciun, but here we do not require that every species suffer a degra- 
dation reaction. This is especially important in studying biochemical 
networks, for which it is rare to have all species degrade. 

1 Introduction 

Although chemical reaction networks can be highly intricate and their in- 
duced differential equations can be extraordinarily complex, it is typically 
the case that dynamical behavior resulting from them is tame. 
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Surprisingly often, for example, there is precisely one stationary state, 
even when kinetic parameters range over all feasible values. This is in con- 
trast to the behavior of complex systems of differential equations in general, 
for which multiple stationary states are common and, in fact, are to be ex- 
pected. Viewed against a background of differential equations generally, then, 
those that derive from chemical reaction networks appear to admit bistable 
behavior far less often than might normally be supposed. This dullness of 
dynamics seems extant across very large classes of reaction networks and per- 
sists robustly, even against large changes in parameter values within kinetic 
models for the rates of individual chemical reactions. 

The implications for the design of biological circuitry are especially com- 
pelling. Evolution has certainly fashioned special reaction network modules 
with dynamical behavior sufficiently rich as to underlie, for example, the 
bistable switching devices or the oscillatory timekeeping mechanisms neces- 
sary for the machinery of life. But these special actors probably require a 
fairly stable and reliable, yet dynamic, chemical background against which 
they can play their roles. In this sense, biology-in-the-large would seem to 
require that only exceptional reaction networks have the capacity to exhibit 
very rich dynamics, while most others execute their functions in steadier, less 
surprising ways. 

To understand, in a general way, the design principles by which biology 
operates, it therefore becomes important to understand the perhaps subtle 
distinctions between reaction network architectures that, on one hand, might 
underlie specialized dynamically-rich devices and, on the other hand, archi- 
tectures that, by their very nature, enforce duller, more restricted behavior 
despite what might be great intricacy in the interplay of many species, even 
independently of values that kinetic parameters might take. 

This article is about a large class of chemical reaction networks, called 
concordant networks, for which a certain dullness of behavior is enforced 
for quite general kinetics. (Concordance is a property of the network itself, 
divorced from any assignment of a kinetics.) Whether a particular network 
is concordant can be readily checked, most easily with help of software of the 
kind made available in [24j, and, as we will show in a separate article, it can 
sometimes be affirmed by inspection of what has come to be known as the 
network's Species- Reaction Graph [ll[T0lll2]. (See Remark [2.11 ) 

Taken with what we call a weakly monotonic kinetics — essentially a 
kinetics in which an increase in the rate of a particular reaction requires that 
there be an increase in the concentration of at least one of its reactant species 



2 



— a concordant network has the following properties (among others): 

(i) Not only can the resulting kinetic system not admit two distinct stoi- 
chiometrically compatible positiv^ equilibria, it cannot even admit two dis- 
tinct stoichiometrically compatible equilibria, at least one of which is positive. 
This is to say that if the resulting kinetic system admits a positive equilib- 
rium, it cannot admit any second equilibrium, positive or otherwise, that is 
stoichiometrically compatible with the first. 

(ii) If the network is weakly reversible and conservative, and if the ki- 
netics is continuous, then among all compositions that are stoichiometrically 
compatible with some positive composition — including "boundary composi- 
tions" (i.e., those for which at least one species concentration is zero) — there 
will indeed be precisely one equilibrium, and it is positive. In fact, the bound- 
ary of the set of all such compositions is repulsive against approach from a 
positive composition in the following sense: At every boundary composition, 
the production rate of each species with concentration zero is non-negative, 
and there is at least one such species whose production rate is strictly pos- 
itive. Moreover, in at least some instances — when the kinetics is smooth 
and the network's "fully open extension" is also concordant (Section [7]) — 
we can assert that the unique equilibrium has a degree of stability: every 
real eigenvalue associated with it is negative. 

In describing these features of concordant networks, we have only told 
part of the story that will unfold later on, but it is perhaps already clear 
that, for concordant networks, the underlying network structure enforces a 
very circumscribed kind of behavior. No matter how intricate a concordant 
network might be — and concordant networks can be highly intricate — it 
cannot, for example, provide the basis for a bistable switch involving two 
equilibria, at least one positive, so long as the kinetics is weakly monotonic. 

This is not to say that concordance in a network precludes all kinds of rich 
behavior over all possible choices of weakly monotonic behavior. It does not. 
(See Remarks 18.51 and 18. 6[ ) What is surprising, however, is that concordance 

— or the lack of it — tells us so much. 

The remainder of this article is organized as follows: In Section [2] we 
discuss some earlier, related results and, in particular, connections to the 
striking work of Banaji and Craciun [21 II]- In Section [3] we provide the 
chemical reaction network theory ideas necessary for the rest of the article. 

"'^A positive composition is one in which all species concentrations are strictly positive. 
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In Section m we provide the definitions of concordance, weakly monotonic 
kinetics, and injectivity of tlie species formation rate function. In the same 
section we show in a very simple way that for a concordant reaction network, 
taken with weakly monotonic kinetics, the species formation rate function 
is injective. A direct consequence of this is the impossibility of two distinct 
stoichiometrically compatible equilibria, at least one of which is positive. In 
fact, we show that, for weakly monotonic kinetics, injectivity is, in a certain 
sense, synonymous with concordance: A reaction network has injectivity in 
all weakly monotonic kinetic systems derived from it if and only if the network 
is concordant. With the results of Section H] behind us, we return briefly in 
Section O to a more fully articulated discussion of connections to other work. 

In Section E] we discuss persistence properties of kinetic systems derived 
from concordant networks, in particular boundary behavior of the species 
formation rate function. Most of the results obtained there rely only on very 
basic, natural features of kinetic rate functions and do not invoke further 
requirements such as weak monotonicity. 

Section [7] is about connections between concordance of a network and 
concordance of a network's "fully open extension" — roughly, the network 
obtained by requiring that every species s of the original network suffers a 
degradation reaction s 0. It is in this section that we begin to make state- 
ments about eigenvalues associated with equilibria. In particular, we show 
how network concordance enforces a certain degree of stability of positive 
equilibria for all (differentiably) monotonic kinetics. By way of contrast, 
in Section |8] we discuss consequences of network discordance beyond those 
stated in Section HI in particular the guarantee of unstable positive equilibria 
for at least certain instances of differentiably monotonic kinetics. Among 
other things we show that for every weakly reversible discordant network 
there is a differentiably monotonic kinetics such that the resulting kinetic 
system admits a positive unstable equilibrium. 

In Section [9] we discuss some variations of the definition of concordance 
and their consequences. There we consider very broad kinetics, in which the 
rate of a particular reaction might even be influenced by the concentration 
of a species that, for the reaction, is neither a reactant nor a product. As an 
important special case, we introduce the notion of strong concordance and 
show that, when the kinetics is similar to the kind admitted by Banaji and 
Craciun [1], then strong concordance suffices for ensuring injectivity of the 
species-formation rate function (and the absence of multiple stoichiometri- 
cally compatible equilibria, at least one of which is positive). This is also 
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discussed in a preliminary way in Section [21 which we turn to next. 

2 Some background and a look ahead 

Chemical reaction network theory has already been able to delineate large 
classes of networks for which dull behavior is ensured, at least with respect 
to certain dynamical features (e.g., the absence of multiple positive stoichio- 
metrically compatible steady states), mostly when the individual reaction 
rates are presumed to be governed by mass action kinetics. In this section 
we discuss two quite separate strands of reaction network theory, and then, 
by way of a preview, we briefly discuss the standing of concordance theory 
in relation to these. 

2.1 Deficiency-oriented theory 

The earliest results were centered about the classification of reaction net- 
works by means of a non-negative integer index called the deficiency. Thus, 
there are networks of deficiency zero, of deficiency one, and so on. The de- 
ficiency is not a measure of a network's size. In fact, a very large complex 
network can have a deficiency of zero. Nevertheless, the Deficiency Zero 
Theorem indicates the sense in which all deficiency zero networks taken with 
mass action kinetics give rise to quite dull dynamics, including the absence 
of multiple stoichiometrically compatible positive equilibria, the absence of 
periodic positive composition trajectories, and the absence of an unstable 
positive equilibrium. The Deficiency One Theorem describes a large class of 
mass action networks for which multiple stoichiometrically compatible pos- 
itive equilibria are precluded, regardless of rate constant values. Still other 
parts of deficiency-oriented theory translate questions about a mass action 
network's capacity for multiple stoichiometrically-compatible positive equi- 
libria — an apparently nonlinear problem — into question about systems of 
linear inequalities. For a survey of some deficiency-oriented reaction network 
theory see [HHIS] . The deficiency-oriented theory continues to evolve . 

2.2 Theory aimed at fully-open reactors 

A rather different strain of reaction network theory began with the Ph.D. 
work of Paul Schlosser [271I2S] and the advent of what was called the Species- 
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Complex- Linkage Graph (SCL Graph). When a reaction network's SCL 
Graph satisfies certain conditions, involving its cycles, then the hypothe- 
sis of a certain technical proposition is also satisfied, and this proposition, in 
turn, ensures that multiple positive steady states are impossible. 

It is important to understand that this work was tailored very specifically 
to determining when a reaction network, taken with mass action kinetics, has 
the capacity to admit multiple positive steady states in the context of what 
chemical engineers call the (isothermal) continuous flow stirred tank reactor 
(GFSTR), in which all species are present in an effluent stream. In reaction 
network theory terms, one studies such reactors by adjoining to the network 
of "true" chemical reactions additional "degradation reactions" of the form 
s 0, one for each species s, to account for the efflux of species s in the 
reactor's effluent [16]. (There are sometimes additional reactions of the form 
s to account for infusion of species s in the feed stream.) Thus, work 
originating in the Schlosser thesis was aimed specifically at understanding the 
capacity for multiple positive equilibria for fully-open reactors, characterized 
by reaction networks in which there is a degradation reaction corresponding 
to each and every species. 

This work was later broadened in [8l[T0l[T2] by Craciun and Feinberg 
to give results which are, in many respects, more powerful, using different 
theoretical underpinnings, but again for mass action kinetics and again for 
reaction networks in which there is a degradation reaction corresponding to 
each species. If the kinetics is mass action, and, for the network at hand, it is 
the case that the Jacobian of the species-formation-rate function is nonsingu- 
lar at every positive composition and for every assignment of rate constants, 
then it was shown that the species-formation-rate function is injective (rela- 
tive to positive compositions) for every assignment of rate constants. In this 
case, multiple positive equilibria become impossible. 

Two techniques were deployed to determine when a network satisfies 
"the Jacobian condition." One was analytical [8], relying on a symbolic 
determinant expansion of the Jacobian with the aim of establishing that, 
for the network under study, every term in the expansion is non-negative 
with at least one positive (whereupon the Jacobian must be nonsingular) . 
The other [T0l[T2], building on the first, was graphical: When the network's 
Species-Reaction Graph (SR Graph) satisfies certain mild conditions, closely 
resembling the earlier ones for the SCL Graph, then the Jacobian condition 
is satisfied, and the species-formation-rate function is injective. 
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Remark 2.1. For the benefit of readers already familiar with results based 
on the Species-Reaction Graph, we reiterate those conditions here: Every 
even cycle of the network's Species- Reaction Graph is an s-cycle and no two 
even cycles have a species-to-reaction intersection. In a separate article we 
intend to show that a "normal" reaction network (Section [7]) is concordant 
— in fact, strongly concordant — if these same conditions are satisfied by 
the network's Species-Reaction Graph. 

In remarkable and highly surprising subsequent work, Banaji and Craciun 
[3111] were able to extend results in (TUlIISj to a far broader class of kinetics, 
hut again for networks in which there is a degradation reaction for every 
species. The kinetics studied by Banaji and Craciun — which they call 
nonautocatalytic (NAC) kinetics — is required to have the property that, for 
a particular reaction, an increase [decrease] in the concentration of one of its 
reactant [product] species — keeping all other concentrations fixed — cannot 
result in a decrease of the reaction rate. These conditions are expressed in 
an essential way in terms of derivatives, so, for them, the kinetics is required 
also to have a degree of smoothness. 

The Banaji- Craciun work built upon earlier work by Banaji and co- 
workers [5], in which it was shown that if the so-called stoichiometric matrix 
of a reaction network has the "SSD property" — i.e, the property that all of 
its square sub-matrices are either singular or else "sign-nonsingular" — and 
if there is a degradation reaction for every species, then, when the kinetics is 
NAC, the species formation rate function is injective (and multiple equilibria 
are excluded). 

Banaji and Craciun jl] then established, surprisingly, that the Species- 
Reaction Graph properties shown to ensure injectivity for certain networks 
in the mass-action case also suffice to ensure that the stoichiometric matrix 
has the "SSD property," which implies in turn that, for fully open reactors, 
one has injectivity of the species-formation rate function even for the very 
broad class of NAC kinetics. 

The requirement that there be a degradation reaction for every species 
figures heavily in the Banaji- Craciun result and in the earlier work by Banaji 
and co-workers: It plays an essential role in connecting the "SSD property" 
of the stoichiometric matrix, via the Gale-Nikaido theorem [20], to injectivity 
of the species-formation rate function. In contrast to the deficiency-oriented 
results described in the preceding subsection, this strong reliance on the 
fully-open setting is a feature common to all of the research described in this 
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subsection, beginning with the Schlosser work, reflecting its origins in the 
study of continuous-flow stirred-tank reactors. 

2.3 A look ahead: stoichiometry, concordance, and 
strong concordance 

In part, fully-open reactors become easier to study because, in that setting, 
stoichiometric constraints become far less of a concern. Reaction networks 
typically studied in biology are not fully open — they rarely contain a "degra- 
dation reaction" for every species — and, for such networks, stoichiometric 
balances come into play more forcefully. Thus, for example, when asking 
about the possibility of multiple equilibria, it is usually sensible to ask only 
about the possibility of two distinct equilibria that are stoichiometrically 
compatible] see Section [31 For this reason, the fully-open results have limited 
range. (Nevertheless, they have been used to infer the absence of multiple 
stoichiometrically compatible equilibria in more general settings [9l[II], but 
the most incisive of the resulting theorems [11] have heretofore been restricted 
to mass- action kinetics.) 

Unlike the more narrowly focused fully-open theory described in §2.2[ the 
deficiency-oriented theory of §2.11 was fashioned to deal with stoichiometric- 
compatibility issues from the outset. For example, when the kinetics is 
mass-action, the Deficiency Zero and Deficiency One Theorems assert the 
impossibility of two distinct stoichiometrically compatible positive equilibria, 
regardless of rate constant values, for any reaction network that satisfies their 
hypotheses. No degradation reactions are required to be present. 

As we shall see in Section HI the definition of network concordance has 
a built-in pre-disposition to stoichiometric compatibility; it is there from the 
very beginning. Only in Section [7] will we be concerned with the concordance 
of a reaction network's "fully open extension" — that is, the network formed 
by adjoining to the original network a full supply of degradation reactions 
— but this is because we can sometimes infer from concordance of the fully 
open extension properties of the original network. (This is similar in spirit 
to [11], but here there is no requirement that the kinetics be mass-action.) 

In a sense, then, the theory of concordant networks shares with the 
deficiency-oriented theory ( §2.ip its indifference to the "fully open require- 
ment," but it also shares with the fully open theory ( §2.2p an affinity for 
special inferences that can be made when there is a full supply of degrada- 
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tion reactions. 

We close this section with an explanation of our choice to defer, com- 
pletely, the notion of strong concordance to the article's end: To the extent 
that they can be compared, the very broad NAC kinetics admitted by Banaji 
and Craciun is significantly broader than the (also broad) weakly monotonic 
kinetics we study in most of this article. (Both subsume mass action kinetics.) 
We show among other things in Section H] that, when the kinetics is weakly 
monotonic, concordance of a reaction network suffices for the impossibility 
of multiple stoichiometrically compatible equilibria. However, in Section [9] 
we show that the same can also be asserted for what we call two-way weakly 
monotonic kinetics, which is very similar NAC kinetics, provided that the 
network is strongly concordant. (In both cases, there is no requirement that 
the reactor be fully open.) 

We chose to center this article around the weaker and broader notion of 
concordance because it already gives rise to a variety of interesting and im- 
portant consequences, some having nothing at all to do with either choice of 
kinetics. By deferring to the article's end the discussion of strong concordance 
and its implications for kinetics of the kind studied by Banaji and Craciun, 
we tried to keep the already powerful consequences of simple concordance 
clearly at center stage. 



3 Reaction network theory preliminaries 

Here we provide the reaction network theory concepts required for stating 
and proving the results in this article. Much of this section is also contained 
in a more extended form, with more discussion and examples, in [TH]. We 
begin with notational conventions. 

3.1 Notation 

In the study of networks, it is sometimes useful to consider vectors whose 
components are real numbers associated with the nodes or with the edges 
of a directed graph. Thus, one typically has a set / (of edges, for example) 
and the components are numbers of the form {xi : i e I}. To vectorize such 
an assignment by working in R^, where N is the number of elements in 
/, requires that one impart an artificial and sometimes awkward cardinal 
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order to /'s elements. Instead, it is useful to simply associate {xi : ^ e 7} in 
the obvious way with a vector x in the vector space of real- valued functions 
having domain /. 

When / is a finite set, we denote the vector space of real- valued functions 
with domain I by W. The subset of W consisting of those functions that 
take only positive (nonnegative) values is denoted R+ (R+). For x & W and 
i e I, the symbol Xi denotes the value assigned to i by x. 

For each x e W , the symbol exp(a;) denotes the element of M+ defined by 

(exp(x))i := expxi, 'ii^I. 

For each x e and for each ^ € M^., the symbol x^ denotes the real number 
defined by 

with the understanding that 0° = 1. For each x, x' e W , the symbol x o x' 
denotes the element of defined by 

(x o x')i ■■- Xix'i, Vi e /. 

For each z e M:(^, the symbol - denotes the element of M.^ defined by 




For each i e /, we denote by oji the element of such that {u;i)j = 1 whenever 
j = i and = whenever j i. 

The standard basis for W is the set {ui e W : i e 1} . Thus, for each x € W , 
we have the representation x = The standard scalar product in W 

is defined as follows: If x and x' are elements of W , then 

tJU JU / lAjq^JUj^m 

i€l 

Note that the standard basis of W is orthonormal with respect to the stan- 
dard scalar product. It will be understood that carries the standard scalar 
product and the norm derived from the standard scalar product. It will also 
be understood that W carries the corresponding norm topology. 

Whenever ^ is a linear subspace of W, we denote by C/-"- the orthogonal 
complement of U in W with respect to the standard scalar product. 
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By the support of x , denoted supp x, we mean the set of indices i € I 
for which Xj is different from zero. When is a real number, the symbol 
sgn (^) denotes the sign of ^. When x is an element of M-^, sgn(a;) denotes 
the function with domain / defined by 

(sgn(x))i := sgnxi, Vie/. 
3.2 Some definitions 

We will use the reaction network displayed in (JT]) to motivate some of our 
definitions. Following Horn and Jackson [22], we call the objects at the 
heads and tails of the reaction arrows — 2A, B,C,C + D and i? in ([1]) — the 
complexes of the network. In this way, we can view the network as a directed 
graph, with complexes playing the role of the vertices and reaction arrows 
playing the role of the edges. 

2A ^ B 

C (1) 

C + D^E 

Remark 3.1. We shall frequently be working in M-^, where is the set of 
species in a network. In this special case, it is advantageous to replace sym- 
bols for the standard basis of R-^' with the names of the species themselves. 
Thus, if the species in the network are given by =5^ = {A, B, C, D, E} then a 
vector such as uc + uJd ^ can instead be written as C + D, and 2ua can 
be written as 2 A. In fact, M.-^ can then be identified with the vector space 
of formal linear combinations of the species. In this way, the complexes of a 
reaction network with species set ^ can be identified with vectors in M-^. 

Definition 3.2. A chemical reaction network consists of three finite sets: 

1. a set ^ of distinct species of the network; 

2. a set ^ c M.-^ of distinct complexes of the network; 

3. a set ^ c ^ X ^ of distinct reactions, with the following properties: 

(a) {y,y)iM for any |/ € "T; 
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(b) for each y ^ ^ there exists y' ^ ^ such that {y,y') e ^ or such 
that {y',y)e^. 

If {y,y') is a member of the reaction set we say that y reacts to y', 
and, following the usual notation in chemistry, we write y ^ y' to indicate 
the reaction whereby complex y reacts to complex y'. We call the complex 
situated at the tail of a reaction arrow the reactant complex of the corre- 
sponding reaction, and the complex situated at the head of a reaction arrow 
the reaction's a product complex. 

The set of species of the network depicted in ([1]) is ^ = {A,B,C,D,E}. 
The set of complexes of the network is = {2 A, B,C,C + D, E}. The set of 
reactions of the network is ^ = {2A B, B ^ C,C ^ 2A,C + D ^ E, E ^ 
C + D). 

The diagram in ([T]) is an example of a standard reaction diagram: each 
complex in the network is displayed precisely once, and each reaction in the 
network is indicated by an arrow in the obvious way. In very few places 
in this article we will refer to the linkage classes of a reaction network, 
which for our purposes can be identified with the connected components 
of the network's standard reaction diagram. Thus, in network ([T]) there are 
two linkage classes, containing, respectively, the complexes {2^4,5,(7} and 
{C + For a formal definition of a linkage class see [I5| . 

More important is the idea of weak reversibility. The following definition 
provides some preparation. 

Definition 3.3. A complex y ^ "io ultimately reacts to a complex e ^ if 
any of the following conditions is satisfied: 

1. y^y'eM- 

2. There is a sequence of complexes y(l),y(2), . . . ,y(k) such that 

y^y{l)^y{2)^...^y{k)^y'. 

In our example, the complex 2A ultimately reacts to the complex C, but 
the complex C does not ultimately react to the complex C + D. 

Definition 3.4. A reaction network {^,^,=^}is called weakly reversible if 
for each y, y' e y' ultimately reacts to y whenever y ultimately reacts to y'. 
A network is called reversible ii y' y e ^ whenever y y' ^ 
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Network ([T]) is an example of a weakly reversible reaction network that 
is not reversible. Note that any reversible network is also weakly reversible. 
Note also that whenever a weakly reversible reaction network is displayed 
as a standard reaction diagram, every arrow in the diagram belongs to a 
directed cycle of arrows. 

Definition 3.5. The reaction vectors for a reaction network l^,*^,^} are 
the members of the set 

{y'-yeR^:y^y'e^}. 

In our example the reaction vector corresponding to the reaction 2A B 
is B - 2A, the reaction vector corresponding to the reaction C + D ^ E is 
E - C - D, and so on. 

Definition 3.6. The stoichiometric suhspace S oia. reaction network ^} 
is the linear subspace of M-^ defined by 

S := span {y' -ye M'^ : y ^ y' e M} . (2) 

We note that, for a reaction network {^,'^,^} , the stoichiometric sub- 
space S will often be a proper subspace of M.-^ . In other words, it will often 
be the case that the dimension of S will be smaller than the number of 
species in the network. For example, in network ([T]) we have dim 5* = 3 while 
dimR^ = = 5. In fact, the stoichiometric subspace will be a proper 

subspace of M'^ whenever the network is conservative: 

Definition 3.7. A reaction network is conservative whenever 

the orthogonal complement 5*-"- of the stoichiometric subspace S contains a 
strictly positive member of M'^: 

S^nRf ^ 0. 

Network ([1]), for example, is conservative: it is easy to verify that the 
strictly positive vector (A + 2B + 2C + D + 3E) e W-f is orthogonal to each of 
the reaction vectors of ([1]). 

If {S^, 3?} is a reaction network, then a mixture state will generally be 

y 

represented by a composition c e , where, for each s e we understand 
Cs to be the molar concentration of s. By a positive composition we mean a 
strictly positive composition — that is, a composition in M.-^ . 
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Definition 3.8. A kinetics ^ for a reaction network is an as- 

signment to each reaction ^/ e ^ of a rate function ^y^y' IR+ 
such that 

J(fy^yi{c) > if and only if suppy c suppc. 

Definition 3.9. A kinetic system {S^ , J(f} is a reaction network ^, £i} 
taken with a kinetics for the network. 

Remark 3.10. In chemical reaction network theory the rate functions are 
generally presumed at the outset to be continuous, but we have not insisted 
on that here. This is because some of the main results in this article (in 
particular those having to do with injectivity) do not require it. For certain 
other formally stated results, we will sometimes say specifically that we are 
assuming continuity or differentiability. On the other hand, when we discuss 
in an informal way the differential equations associated with a kinetic sys- 
tem, it will be understood implicitly that there is smoothness sufficient to 
ensure the that the commonly presumed properties of differential equations 
are present. 

In this article we shall often refer to mass action kinetics and to mass 
action kinetic systems. Both are formally defined below: 

Definition 3.11. A kinetics J(f for a reaction network {=5^,^,^} is mass 
action if, for each reaction y ^ y' e M, there is a positive number ky^y' such 
that the rate function J^^y' takes the form 

J(^^y'{c) = ky^ylCP . (3) 

The positive number ky^yi is the rate constant for reaction y ^ y' . 

Definition 3.12. A mass action kinetic system is a reaction network taken 
together with a mass action kinetics for the network. 

Definition 3.13. The species formation rate function for a kinetic system 

y 

{S^ , with stoichiometric subspace 5" is the map / : S defined 

by 

/(c) = X,^y'(c)(y'-y). (4) 
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For a kinetic system {^,^,^, J^} whose underlying network is our ex- 
ample network ([T]), the species formation rate function has the following 
species-wise form: 



/a(c) = -2je2A^B{c) + 2jeB^2A(c) + 2jec^2A{c), (5) 

/c(c) = -^^2a(c) + J^^c(c) - ^c+d^e{c) + J(^e^c+d{c), 

fsic) = -J(^E^C+d{c) + J(fc+D^E{c). 

Definition 3.14. The differential equation for a kinetic system with species 
formation rate function /(•) is given by 



c = /(c). (6) 

Let ^} be a kinetic system. From equations ([2D, Q, and (jS]) 

we observe that the vector c will invariably lie in the stoichiometric subspace 
S of the network . Thus, the difference of any two compositions 

y y 

c e and c' e IR+ that lie along the same solution of (E]) will always reside 
in 5*. This motivates the following definition: 

Definition 3.15. Let be a reaction network with stoichiometric 

y 

subspace S. Two compositions c and c' in are called stoichiometrically 
compatible if c' - c € 5". 

We note that stoichiometric compatibility is an equivalence relation. As 

— y 

such, it partitions R_,_ into equivalence classes that we call stoichiometric 
compatibility classes. Thus, the stoichiometric compatibility class containing 

y 

an arbitrary composition c, denoted (c + S*) n ]R_|_ , is given by 

(c + ^)nlf = {c'€lf :c'-c€>9}. (7) 
We observe, as the notation suggests, that {c + S) n M'f is the intersection of 

y 

]R_i_ with the parallel of S containing c. 

A stoichiometric compatibility class will typically contain a wealth of 
(strictly) positive compositions. We say that a stoichiometric compatibility 
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class is nontrivial if it contains a member of ]R':f . To see that a stoichiomet- 
ric compatibihty class can be trivial, consider the simple reaction network 
A + B C, and let c be the composition defined by ca = 1,cb = 0,cc = 0. 
Then the stoichiometric compatibility class containing c has c as its only 
member. 

Definition 3.16. An equilibrium of a kinetic system is a 
y 

composition c e for which /(c) = 0. A positive equilibrium of a kinetic 
system {^,^,,^, J^} is an equilibrium that lies in Rf. 



In light of Definition [3]8l a kinetic system can admit a positive equilibrium 
only if its reaction vectors are positively dependent: 

Definition 3.17. The reaction vectors for a reaction network {J^ ,M} 
are positively dependent if for each reaction y ^ y' e M there exists a positive 
number Ky^y' such that 

E S-y'(2/'-2/) = 0. (8) 

Remark 3.18. For any weakly reversible network, the reaction vectors are 
positively dependent [T5] . 



Occasionally we will want to consider changes in the value of the species 
formation rate function corresponding to small departures from a positive 
composition c* toward nearby compositions that are stoichiometrically com- 
patible with c*. Thus, for a kinetic system {=5^,^,^,J^} with stoichio- 
metric subspace S, with smooth reaction rate functions, and with species 

— .y 

formation rate function / : ^ S", we will want to work with the derivative 
df{c*) S S, given by 



de 



, ^aeS. (9) 

=0 

In this case, we say that c* € is a degenerate equilibrium if c* is an 
equilibrium and if, moreover, df(c*) is singular. 
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4 Concordant networks, weakly monotonia 
kinetics, and injective kinetic systems 

In preparation for the definition of concordance we consider a reaction net- 
work ^} with stoichiometric subspace S c M-^, and we let L : -> S 
be the hnear map defined by 

La= Oiy^y'iy' -y)- (10) 

Definition 4.1. The reaction network {S^ ,'10 is concordant if there do 
not exist an a e ker L and a nonzero a e S having the following properties: 

(i) For each y y' e 3^ such that ay^y' ^ 0, supp?/ contains a species s for 
which sgn CTs = sgna^^j,/. 

(ii) For each y y' e M such that cty^y^ = 0, either as = for all s € supp?/ 
or else suppy contains species sand s' for which sgnas = - sgn (7^', both not 
zero. 

A network that is not concordant is discordant. 

Remark 4.2. Concordance is a network property that is discernible by sign- 
checking in kerL and S. A Windows-based computer program that imple- 
ments such a concordance test is freely available [21]; it is fast for networks 
of moderate size. 

Examples. While network ([1]) is a fairly simple "toy" that is concordant, the 
network depicted in (ITT]) serves as a more intricate and more biologically- 
motivated concordant example. In Reidl et al. [26], network ( !TT|1 served to 
model calcium dynamics in olfactory cilia. The concordance of (fTT]) was 
ascertained by means of the sign-checking procedure implemented in [24] . 

C + AD^E 
B + E^ F^A + E (11) 
B^ D + B D^O 

A biologically motivated example of a discordant network is shown in 
f lT2|) . This network underlies a model of the bacterial two-component sig- 
naling system EnvZ-OmpR [2UH^ . which regulates the number and size of 
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membrane pores in response to changes in osmolarity. 



C + D^E^A + F 
B+F^G^B+D 



(12) 



Definition 4.3. A kinetic system {y ^M^J^) is injective if, for each pair 

y y 

of distinct stoichiometrically compatible compositions c* € and c** e M+ , 
at least one of which is positive, 



Remark 4.4. Clearly, an injective kinetic system cannot admit two distinct 
stoichiometrically compatible equilibria, at least one of which is positive. 

Definition 4.5. A kinetics for reaction network {S^ is weakly 
monotonic if, for each pair of compositions c* and c**, the following impli- 
cations hold for each reaction y y' e ^ such that suppy c suppc* and 
suppy c suppc**: 

(i) J^y^y>{c**) > J^y^yi{c*) => thcrc Is a species s e suppt/ with c** > c*. 

(ii) J^y^yi{c**) - J(fy^yi{c*) c* * = c* for all s e supp?/ or else there 
are species s, s' € suppy with c** > c* and c** < c*,. 

We say that the kinetic system {S^ ,M,J(f} is weakly monotonic when its 
kinetics is weakly monotonic. 

Remark 4.6. Note that if a kinetic system is weakly monotonic and if c** 
is a positive composition, then the implications (i) and (ii) will hold for all 
reactions in the network, even when c* is not strictly positive. 

Remark 4.7. Clearly, every mass action kinetic system is weakly monotonic. 
Note that the definition of weak monotonicity makes no assumptions about 
the continuity or the smoothness of the kinetic rate functions. When, for a 
weakly monotonic kinetics, a reaction-rate function J^^j^'(-) is differentiable, 
there is no requirement that its derivative with respect to c^, s e suppy, be 
strictly positive at all positive compositions. 



^ jey^y,{cn{y'-y)+ E ^y^yimy' -y). 



(13) 
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Proposition 4.8. A weakly monotonic kinetic system {^,^,^,J^} is in- 
jective whenever its underlying reaction network is concordant. In 

particular, if the underlying reaction network is concordant, then the kinetic 
system cannot admit two distinct stoichiometrically compatible equilibria, at 
least one of which is positive. 

Proof. Suppose on the contrary that, for the weakly monotonic kinetic sys- 
tems ,M,J(^}, there are distinct stoichiometrically compatible compo- 
sitions c* and €**, with c** e Rf"', such that 

^ X,^Acniy'-y)= E Xj^y'ic*)iy'-y). (14) 

Clearly, then, a € defined by 

ay^y. := Xj^y'iC**) - Xj^yiC* ) , V ^ ' € ^ (l5) 

is a member of ker L. Because c* and c** are distinct and stoichiometrically 
compatible, the vector a := c** - c* is a nonzero member of S. Note that if, 
for a particular reaction y y', we have o^y^y' > 0, then weak monotonicity 
and Remark 14.61 require that, for some species s e suppy, as ■= c** - c* > 0. 
A similar argument can be mounted when we have cty^y' < 0; in this case, 
there is s € supp y such that ag c** - c* < 0. Finally, if we have cty^y^ - 0, 
then weak monotonicity requires that the numbers {cr^jsesuppi/ are either all 
zero or else there are two of opposite nonzero sign. The existence of a and 
a, so constructed, contradicts the supposition that the network ,M} 
is concordant. □ 

Example. Proposition 14.81 tell us that no weakly monotonic kinetics for the 
complex but concordant network ffTTj) can give rise to multiple stoichiometri- 
cally compatible equilibria, at least one of which is positive. This is consistent 
with the findings in [26] . 

The simple Proposition 14.81 — one of the most important of this article 
— tells us that, for any weakly monotonic kinetic system, injectivity (and 
therefore the absence of distinct stoichiometrically compatible equilibria, at 
least one of which is positive) can be precluded merely by establishing con- 
cordance of the underlying reaction network, perhaps with the assistance of 
a concordance-checking computer program such as the one made available 
in 121]. As we shall see in a subsequent article, the Species- Reaction Graph 
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can sometimes also serve to establish network concordance. (See Remark 

[231) 

The converse of Proposition I4.8l is not true. That is, there can be a weakly 
monotonic kinetic system — in fact, a mass action system — that is injective 
even when its underlying reaction network is not concordant. On the other 
hand, the following proposition resembles a converse; proof by construction 
is provided in Appendix A. 

Proposition 4.9. Let {^,^,^}6e a discordant reaction network. Then 
there exists for it a weakly monotonic kinetics such that the resulting 
kinetic system {^,^,,^, J^} is not injective. 

Remark 4. 10 (Example) . Proposition 14.91 asserts the existence of a weakly 
monotonic kinetics for the discordant network (fT2|) that gives rise to a non- 
injective kinetic system. In fact, it can be shown by direct computation that 
any mass-action system based on ( IT^ possesses at least two stoichiometically 
compatible equilibria [21] — one positive and one non-positive — thereby 
violating injectivity. 

Taken together. Propositions 14. 81 and 14.91 can be summarized in the follow- 
ing theorem — a theorem that tells us that the class of concordant reaction 
networks is precisely the class of networks that are injective for every assign- 
ment of a weakly monotonic kinetics: 

Theorem 4.11. A reaction network has injectivity in all weakly monotonic 
kinetic systems derived from it if and only if the network is concordant. 

We discuss more subtle connections between concordance and injectivity 
in the next section. 

5 A remark about connections to earlier work 

In [Hllin] and [12] attention was focused on mass action kinetic systems in 
which the underlying reaction network contains a "degradation reaction" of 
the form s -> for each species s in the network. (Here such networks 
are called fully open.) In this context it was shown that there are certain 
fully open networks having the property that all mass action kinetic systems 
deriving from them are injective — that is, regardless of rate constant values. 
These were called injective (mass action) networks in [5|[TU] . and it was shown 
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that a network is injective in this (mass action) sense if its Species- Reaction 
Graph satisfies the mild conditions stated in Remark 12.11 

With this in mind, it is not unreasonable to suppose that a variant of 
Proposition 14.91 might be true: A network can have the property that every 
mass action system deriving from it is injective, regardless of rate constants, 
only if the network is concordant. 

In fact, this is false. We show in f lTB]l an example of a fully open reaction 
network that is not concordant but for which any mass action system derived 
from it is injective. Injectivity follows from a theorem in [7ji that invokes 
Species- Reaction Graph conditions different from those in Remark 12.11 On 
the other hand, it can be established that the network is not concordant by 
a choice of a and a that satisfy the conditions in Definition 14.11 

A + B^C 2A + B^D A^O^B (16) 

C D 

The example tells us that "network injectivity" in the mass action case 
should not be conflated with network concordance. In particular, the exam- 
ple tells us that there are fully open networks for which particular Species- 
Reaction Graph conditions, different from those in Remark \2. 1\ might ensure 
^^mass action network injectivity" but for which network concordance need 
not obtain. 

In a separate article, however, we will show that when a fully open net- 
work's Species-Reaction Graph satisfies the specific conditions stated in Re- 
mark 12.11 — conditions that ensured injectivity in the mass action setting 
considered in ^[101 [12] and in the far more general NAG kinetics setting 
considered by Banaji and Graciun 0,0] — then those same Species-Reaction 
Graph conditions will ensure that the network itself, divorced from any kinet- 
ics at all, is concordant and, in fact, strongly concordant ( Section [^. Thus, 
when the kinetics is weakly monotonic or two-way weakly monotonic (Section 
E]), injectivity and the absence of multiple equilibria are ensured. 

In the fully open network setting, then, the Species-Reaction Graph con- 
ditions stated in Remark |2. II will imply injectivity for large classes of kinetics 

similar theorem appeared earlier in [27l[25. In the case of ([T6)) . injectivity follows 
from the fact that the Species-Reaction Graph has just one cycle, and it is a c-cycle — 
that is, a cycle whose edge set is the union of c-pairs. 
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through the very simple Propositions 14.81 and 19.61 While earlier work based 
on these conditions appealed to Jacobians, special determinant attributes, and 
the Gale-Nikaido theorem to infer injectivity, the elementary Propositions \4 ■ 8\ 
and \9.6\ make those appeals unnecessary. Moreover, we shall see in Section [7] 
how properties of the Species-Reaction Graph will often imply concordance 
of a network (and, therefore, injectivity of any weakly monotonic kinetic 
system derived from it) even when the network is not fully open. 

6 Persistence in weakly reversible concordant 
kinetic systems: Boundary behavior and 
the existence of positive equilibria 

It is our aim in this section to establish two related properties inherent to 
kinetic systems, not necessarily weakly monotonic, that derive from weakly 
reversible concordant reaction networks. 

These properties are connected to "persistence" questions: Viewed from 
an ecological perspective, one might be interested to know for a kinetic sys- 
tem whether there will be a positive equilibrium — that is, an equilibrium 
in which all species coexist. And one might also want to know whether a 
population that begins at a state in which all species are present can evolve 
to a different one — in particular to an equilibrium — in which one or more 
species are extinct (i.e., absent). For kinetic systems in which the underlying 
network is weakly reversible, questions like these become especially challeng- 
ing. 

As we shall see, however, a weakly reversible concordant reaction network 
invariably has pleasant properties that make these questions far easier to 
resolve: 

First, we examine behavior of the species formation rate function on the 
boundary of a nontrivial stoichiometric compatibility class. (Recall that by 
a nontrivial stoichiometric compatibility class we mean one that contains a 
strictly positive composition.) In particular, we show that, on the boundary, 
there can be no equilibria at all. In fact, at every boundary composition c, the 
species formation rate vector /(c) points into the stoichiometric compatibility 
class in the following sense: There is at least one species s* for which Cg* = 
while Cs* = fs*{c) > 0. (Thus, the inward-pointing boundary vector field 
thwarts an approach to the boundary along a trajectory that originates at a 
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strictly positive composition.) 

Second, we show that if the network is conservative and the kinetic rate 
functions are continuous then each nontrivial stoichiometric compatibihty 
class contains at least one positive equilibrium. If, in addition, the kinetics is 
weakly monotonic, then there is precisely one equilibrium in each nontrivial 
stoichiometric compatibility class, and it is positive. (Of course, in the weakly 
monotonic case, uniqueness already follows from Proposition 14.81 ) 

6.1 Preliminaries: reaction transitive compositions 



We begin with some old reaction network theory ideas, derived from 



that are independent of concordance. The following definition is intended to 
convey a simple idea: A species set S^* c ^ is reaction-transitive if, whenever 
all the reactant species of a reaction are present in then so too are that 
reaction's product species. 

Definition 6.1. For a reaction network , a species set ^* c ^ 

is reaction-transitive if, for all y ^ y' e 



We say that a composition c € is reaction-transitive if supp c is reaction 
transitive. 

Clearly every positive composition is reaction-transitive. It follows from 
properties presumed of a kinetics in Section [3] that if c is a composition on 

the boundary of and if species s is "missing" at composition c (i.e., if 
Cs = 0), then the formation rate of s at composition c is not negative (i.e., 
/s(c) > 0). The following lemma tells us that if c is not reaction-transitive, 
then there must be a missing species which, at composition c, is produced at 
a strictly positive rate. (See [U] or Appendix I in 

Lemma 6.2. Consider a kinetic system with species forma- 

y 

tion rate function /(•). If c is a composition on the boundary o/ ]R_,_ , then 

^See in particular Propositions 5.3.1-5.3.2, Remark 6.1.E, and Appendices I and II 
in [16] . The idea behind Del. 16.11 was important not only in [16] but also in more recent 
papers, where the formulation is phrased differently. A species set is called semi-locking 
in [T] or a siphon in if its complement in the full set of species is reaction-transitive. 




suppyc^* =^ suppy'c^*. 



(17) 
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/s(c) > for each species s such that Cg = 0. Moreover, if c is not reaction- 
transitive, then there is a species s* such that Cg* = and /s*(c) > 0. Thus, 
every equilibrium of the kinetic system is reaction transitive. 

6.2 Boundary behavior for weakly reversible concor- 
dant networks 

With this as background, we are in a position to state a theorem about 
weakly reversible concordant reaction networks. Note that the theorem makes 
mention only of network properties; there is no mention of kinetics. 

Theorem 6.3. For a weakly reversible concordant reaction network, no non- 
trivial stoichiometric compatibility class can have on its boundary a compo- 
sition that is reaction-transitive. 

Before proving Theorem 16.31 we observe that the theorem and Lemma 
16.21 give us the following corollary: 

Corollary 6.4. For a kinetic system that derives from a weakly reversible 
concordant reaction network, no composition on the boundary of a nontriv- 
ial stoichiometric compatibility class is an equilibrium. In fact, at each such 
boundary composition the species formation rate vector points into the sto- 
ichiometric compatibility class in the sense that there is an absent species 
produced at strictly positive rate. 

As might be expected, a species- formation-rate function that, in the sense 
indicated, points inward on the boundary of a non-trivial stoichiometric com- 
patibility class repels an approach to the boundary by a bounded trajectory 
that begins at a strictly positive composition. For a formal proof that its 
cj-limit set is, in fact, disjoint from the boundary when the kinetics is suffi- 
ciently smooth see [2]. 

Remark 6.5. (About the proof of Theorem 1 6. 31) Although Theorem 16.31 is 
about network properties alone, divorced from any mention of kinetics, in its 
proof we will endow the reversible concordant network of the theorem with 
a mass action kinetics in order to produce a contradiction. To begin, we 
will presume that there can be a nontrivial stoichiometric compatibility class 
that has a reaction-transitive composition c on its boundary. Then we will 
endow the network with a mass action kinetics, with rate constants chosen 
to make the boundary composition c an equilibrium. 
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Now suppose we could take for granted that, for every weakly reversible 
network and for every assignment of rate constants, the corresponding mass 
action species-formation rate function admits a positive equilibrium in each 
nontrivial stoichiometric compatibility class. In that case, we would have 
a contradiction of Proposition 14.81 In particular, in the nontrivial stoichio- 
metric compatibility class containing c we would have, for the mass action 
system constructed as described, a positive equilibrium that is stoichiomet- 
rically compatible with the boundary equilibrium c. 

In fact, an unpublished manuscript by Deng et al. [13] does indeed provide 
an argument for the existence, in each nontrivial stoichiometric compatibility 
class, of a positive equilibrium for any mass action system in which the un- 
derlying network is weakly reversible. In the absence of a published version, 
however, we can instead mount a similar but somewhat more complicated 
argument that relies on the same strategy, this time involving construction of 
mass action system that is complex balanced at the boundary composition c. 
This would permit invocation of known properties of complex balanced mass 
action systems, in particular the distribution of their positive equilibria. 

Because the argument that invokes \13l is more straightforward, we give it 
here in the main text. The alternative argument, which does not rely on [13], 
is provided in Appendix B. 

Remark 6.6. In preparation for the nearby proof of Theorem 16.31 (and of the 
proof in the Appendix), it will be helpful if we make the following observation 
[TB] : Consider a weakly reversible network , and suppose that a 

composition c* is reaction-transitive. It is not difficult to see that, if and 
y"^^ are complexes belonging to the same linkage class, then there are two 
possibilities: Either supp y^ and supp y^^ are both contained in supp c* or else 
neither is. 

This is to say that the linkage classes are of two distinct kinds: Relative 
to the reaction-transitive composition c*, a supported linkage classes is one in 
which the support of every complex is contained in suppc*, while an unsup- 
ported linkage class is one in which the support of no complex is contained 
in supp c*. 

If the network is given a kinetics, it follows from properties of a kinetics 
given in Section [3] that, at the reaction-transitive composition c* , all reactions 
associated with a supported linkage class will proceed at strictly positive rate, 
while all reactions associated with an unsupported linkage class will have zero 
rate. In rougher terms, at composition c*, all reactions associated with a 
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supported linkage class will be "turned on" while those associated with an 
unsupported linkage class will be "turned off." Moreover, the "turned on" 
reactions (if there are any) constitute a weakly reversible subnetwork of the 
original network. 

Proof of Theorem \6.3[ Suppose that is a concordant weakly re- 

— y 

versible reaction network and that, contrary to what is to be proved, c* e 
is a reaction-transitive composition on the boundary of a nontrivial stoichio- 
metric compatibility class. Let 

^* :={?/€ ^ : supp y c supp c*} and ^* ■■= {y ^ y' e ^ : supp y c supp c*}. 

Provided that ^* is not empty, {S^,'t^*,^*} is a weakly reversible sub- 
network of the original network. In this case, there are positive numbers 
{Ky^y'}y^y'e.'ig* that satlsfy [in] 

E S-/(2/'-?/) = 0. (18) 

Now choose rate constants {ky^yi}y^yi^_>^* for reactions in to satisfy 

S-^/' = W(c*)^ ^y^y'e^*. (19) 

In any case, for the remaining reactions choose rate constants {ky^yf}y^yf^^sSi* 
to be any positive numbers. Note that for any y y' e ^ \ ^* we will have 
ky^yr(c*)y = because suppy / suppc*. From this, ([TS]) . and (fT^ it follows 
that 

E ky^y,{c*)y{y'-y) + ^ ky^y,{cy{y'-y) = 0. (20) 

This is to say that c* is a (boundary) equilibrium of the mass action sys- 
tem {y,'T^,^,k} so constructed. Because {^,^,^}is weakly reversible, 
[13] asserts, for the same mass action system, the existence of a strictly 
positive equilibrium c** in the nontrivial stoichiometric compatibility class 
containing c*. Because the network is concordant, this contradicts the con- 
clusion of Proposition 14.81 □ 
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6.3 Conservative weakly reversible concordant networks 
and the existence of positive equilibria 

When, in the following theorem, we say that a kinetics is continuous we mean 
that its various reaction rate functions are all continuous. 

Theorem 6.7. // ^ is a continuous kinetics for a conservative reaction 
network {S^ , then the species formation rate function for the kinetic 
system J^} has an equilibrium within each stoichiometric compat- 

ibility class. If the network is weakly reversible and concordant, then within 
each nontrivial stoichiometric compatibility class there is a positive equilib- 
rium. If, in addition, the kinetics is weakly monotonic, then that positive 
equilibrium is the only equilibrium in the stoichiometric compatibility class 
containing it. 

We begin with a lemma that is required for application of a fixed point 
theorem due to Browder [H] . For any kinetic system, the lemma merely makes 
more firmly geometrical the idea that, on the boundary of a stoichiometric 
compatibility class, the species formation rate vector does not point outward. 

Lemma 6.8. Consider a kinetic system with stoichiomet- 

— y 

ric subspace S and species formation rate function / : S. If c is a 

y 

composition on the boundary of R_|_ , then there is a composition c^ in the 

stoichiometric compatibility class (c + S) nR_^ and a positive number A such 
that f{c) = X{ct-c). 

Proof. Choose 9 > sufficiently small as to ensure that 

Cs + 9fs{c) > 0, Vs e supp c. (21) 
For s i suppc we have = and, by Lemma [621 /s(c) > 0, so 

Cs + 9fsic) > 0, ysi supp c. (22) 

Now let c^ ■= c + 9f{c). From ( 12T|) and ( l22l) it follows that is a member 

of M_,_ . Because /(c) is a member of S, it follows that c'f is a member of 
the stoichiometric compatibility class containing c. Taking A = |, we get the 
desired result. □ 
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Proof of Theorem \6. ?[ Let c be a stoichiometric compatibility class. 

y 

That is, there is a composition c# € such that 

K={c* + S)n^^. 

Because the network is conservative, K is compact [22], and 

it is easily seen to be convex. Let g ■■ K ^ (c# + S) be defined by g{c) ■■- 
f{c) + c, where /(•) is the species formation function for the kinetic system 
From Lemma [6.81 it is apparent that for each composition c 
on the boundary of K there is a number A and a composition e K (both 
depending perhaps on c) such that g{c) - c = A(c^ -c). In this case, it follows 
from a theorem of Browder [6] that g{-) has a fixed point in K — that is, a 
composition c* e K such that g{c*) = c* . But then we must have /(c*) = 0, 
which is to say that there is an equilibrium in K. 

If the network is weakly reversible and concordant and if the stoichiomet- 
ric compatibility class K in nontrivial, it follows from Corollary 16.41 that c* 
is not on the boundary of K\ it must be positive. If, in addition, the kinetics 
is weakly monotonic, the fact that there cannot be another equilibrium in K 
follows from Proposition 14.81 □ 



7 Concordance and properties of a network's 
fully open extension; eigenvalues 

By the fully open extension of reaction network {^,^,=^}we mean the 
network obtained from the original network by adjoining to it the species 
degradation reactions {s ^ : s e y} if those reactions are not already 
present. More precisely, the fully open extension of is the net- 

work {^,'^,^1, wherfl 

<^:='^u^u{0} and M := M u {s ^ : s ^ y} . (23) 

In this section we examine implications of concordance in a reaction net- 
work and in its fully open extension. We begin by noting that a network 
need not be concordant even when its fully open extension is concordant. 
The simple network B -^^ A ^ C provides an example. 

^Note that in the definition of ^ we are viewing as a set in M'^ — that is, as a 
surrogate for the standard basis of M-^. See Section E] 
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Remark 7.1 {Why study a network's fully open extension?). There are two 
reasons to examine the concordance of a network's fully open extension. The 
first of these is that concordance of the fully open extension gives partial 
information about the stability properties of equilibria for kinetic systems 
derived from the original network, at least when the kinetics satisfies certain 
modest requirements. This we shall see in Theorem 17.61 and in Section [HI 

The second reason is connected to inferences one can make from the 
original network's Species-Reaction Graph. As we shall see in a subsequent 
article, properties of a network's Species-Reaction Graph will often indicate 
that the fully open extension of a network is concordant. For this reason, we 
will want to know when fully-open-network concordance implies concordance 
of the original network. This is the subject of Theorem 17.41 which also has 
some stability implications, discussed in Section [81 

In a mass action context, connections were made in [TT] between injec- 
tivity of a reaction network and of its fully open extension. Here we will 
be able to generalize those results substantially. By way of preparation, we 
record the following definition: 

Definition 7.2. Consider a reaction network {=5^,"^,^} with stoichiometric 
subspace S. The network is normal if there are q e and 77 € R-* such that 
the linear transformation T : S ^ S defined by 

T^--= E Vy^y'iy*^)iy'-y) (24) 

y^y'e.^ 

is nonsingular, where is the scalar product in M'^ defined by 

X * x' ■■= ^ QsXsX'g. (25) 

Remark 7.3. The normality condition is not difficult to satisfy. In fact, it 
was shown in [11] that every weakly reversible network is normal, regardless 
of the complexes. Weaker sufficient conditions for network normality are also 
given in [TT] . 

We are now in a position to state a theorem that has considerable signifi- 
cance when placed alongside another one, to appear in a forthcoming article, 
that indicates how properties of a network's Species-Reaction Graph can en- 
sure concordance of the network's fully open extension. In that theorem. 
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however, concordance of the original network is left undecided. The follow- 
ing theorem completes the picture when the original network is normal. As 
we shall see, Theorem 17.41 also has some stability consequences. A proof is 
provided in Appendix C. 

Theorem 7.4. A normal network is concordant if its fully open extension 
is concordant. In particular, a weakly reversible network is concordant if its 
fully open extension is concordant. 

So far in this section all considerations were about the network alone, di- 
vorced from a kinetics. Now we will suppose that a network ^} is en- 
dowed with a kinetics J^such that the resulting kinetic system {y, 
admits a positive equilibrium c*. We ask what network concordance — of 
either the network or of its fully open extension — might tell us about eigen- 
values of df{c*) S ^ S, where S is the stoichiometric subspace for the 
network, /(•) is the species formation rate function for the kinetic system, 
and df(c*) is the derivative of /(■) at c*. 

In this case, we will require that the kinetics be "differentiably monotonic" 
at c*. (It is easy to see that every mass action kinetics for a network is 
differentiably monotonic at every positive composition.) 

Definition 7.5. A kinetics Jt^ for a reaction network {S^ is differ- 
entiably monotonic at c* € Rf' if, for every reaction y y' e J^^y'{-) is 
different iable at c* and, moreover, for each species s e 

^^{c*)>0, (26) 

with inequality holding if and only if s e supp y. A differentiably monotonic 
kinetics is one that is differentiably monotonic at every positive composition. 

Theorem 7.6. Let be a kinetic system with stoichiometric 

y 

subspace S and species formation rate function / : M_|_ S. Moreover, sup- 
pose that the kinetics is differentiably monotonic at c* e Mf" . If the underlying 
network {<y,^,^}zs concordant then the derivative df{c*) ■ S ^ S is non- 
singular (whereupon is not one of its eigenvalues). If the network's fully 
open extension is concordant then no real eigenvalue of df(c*) is positive. 

Proof of Theorem 17.61 is given in Appendix [Dl 
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Corollary 7.7. Let {=5^,^,^, J^} be a kinetic system with stoichiometric 

y 

suhspace S and species formation rate function / : ^ S. Moreover, sup- 
pose that the kinetics is differentiahly monotonic at c* e Mf^. If the underlying 
network {^,^,,^} is normal — in particular, if it is weakly reversible — and 
if the network's fully open extension is concordant, then every real eigenvalue 
of df{c*) is negative. 

Proof. If the underlying network is normal then, by virtue of Theorem I7.4[ 
its concordance follows from concordance of the fully open extension. In this 
case Theorem 17. 6l ensures that df{c*) is nonsingular, so cannot be among its 
eigenvalues. By the same theorem, concordance of the fully open extension 
then ensures that all real eigenvalues are negative. □ 

8 Consequences of discordance 

Proposition 14.81 told us that a concordant network, taken with any weakly 
monotonic kinetics, invariably gives rise to an injective kinetic system. In 
this sense, concordance is a network attribute that enforces injectivity for 
all weakly monotonic kinetics. Proposition 14.91 went even further: It told us 
that no discordant network has this same property. 

The observations in this section are in the same spirit, this time as coun- 
terpoints to Theorem 17.61 and its corollary. Those told us that a normal 
network with concordant fully open extension, taken with any differentiahly 
monotonic kinetics, will invariably have some degree of stability at any pos- 
itive equilibrium it might admit, in the sense that every real eigenvalue is 
negative. Here we shall see that no network with discordant fully open ex- 
tension — in particular, no weakly reversible discordant network — can have 
this same property (except in the trivial case for which there are no positive 
equilibria at all). 

Because our interest will be in positive equilibria, our focus will be on 
networks with positively dependent reaction vectors. (Recall from Section [3] 
that a kinetic system can admit a positive equilibrium only if the reaction 
vectors for the underlying network are positively dependent.) Proof of the 
following theorem is provided in Appendix [Dl 

Theorem 8.1. Consider a reaction network with positively dependent reac- 
tion vectors. If the network is discordant, then there exists for it a differ- 
entiably monotonic kinetics such that the resulting kinetic system admits a 
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positive degenerate equilibrium. If the network's fully open extension is dis- 
cordant then there exists for the original network a difjerentiably monotonia 
kinetics such that the resulting kinetic system admits an unstable positive 
equilibrium — in fact, a positive equilibrium associated with a positive real 
eigenvalue. 

Theorem 18.11 has a number of corollaries. In each corollary, it is left 
implicit that the unstable positive equilibrium mentioned can be taken to 
have associated with it a positive real eigenvalue. 

Corollary 8.2. For every discordant normal network with positively depen- 
dent reaction vectors there exists a difjerentiably monotonia kinetics such that 
the resulting kinetic system admits a positive unstable equilibrium. 

Proof. It follows from Theorem 17. 41 that every normal discordant network has 
a discordant fully open extension. Corollary 18.21 then follows from Theorem 
11 □ 



Corollary 8.3. For every weakly reversible discordant network there is a 
differentiably monotonia kinetics such that the resulting kinetic system admits 
a positive unstable equilibrium. 

Proof. A weakly reversible network is normal and, moreover, has positively 
dependent reaction vectors. Corollary 18.31 is then a consequence of Corollary 
[Q □ 

In preparation for the next corollary we point out that a concordant 
weakly reversible network (in fact, a reversible concordant network) can have 
a discordant fully open extension. Network (|27|) provides an example, as 
indicated by computation ^]. (This is a variant of a network considered for 
different reasons in j28].) 

A + B^E C + D^F (27) 
C ^ 2A ^ D ^ B 

The following is a consequence of Proposition 14.81 Corollary 16. 4[ and 
Theorem 18. 1[ 



Corollary 8.4. For every weakly reversible concordant reaction network with 
discordant fully open extension there exists for the original network a differ- 
entiably monotonia kinetics such that the resulting kinetic system (for the 
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original network) has the following properties: There is an equilibrium that 
is positive, unique within its stoichiometric compatibility class, and unstable. 
Moreover, at each point on the boundary of that stoichiometric compatibility 
class the species-formation-rate vector points inward in the sense that there 
is an absent species produced at strictly positive rate. 

Remark 8.5. Although network concordance, through Proposition 14. 8[ en- 
sures the uniqueness of a (positive) equilibrium within its stoichiometric 
compatibility class, network (!27|) and Corollary 18.41 tell us that concordance 
by itself cannot also ensure that the sole equilibrium is stable. 

We will show in another article that when a normal (in particular, a 
weakly reversible) network's Species-Reaction Graph satisfies the conditions 
mentioned in Remark 12. ![ concordance of both the network and its fully 
open extension follow. In this case, the hypothesis of Corollary 18.41 cannot 
be satisfied. 

Nevertheless, even when both a network and its fully open extension are 
concordant, and even when the kinetics is differentiably monotonic, we still 
cannot preclude instability of the unique positive equilibrium: Although its 
associated real eigenvalues must be negative ( Corollary 17. 7p . there might also 
be complex eigenvalues with positive real part. Such eigenvalues were found 
to be extant in Reidl et al. [2B], where network ( II ip . taken with a weakly 
monotonic (almost mass action) kinetics, served as a model for calcium os- 
cillations in olfactory cilia|§ As indicated by computation [24J, both network 
ffTTj) and its fully open extension are concordant. 

Thus, it appears that concordance of a network (and of its fully open 
extension) enforce a degree of dull behavior, at least to the extent that they 
preclude multiple stoichiometrically compatible equilibria (and positive real 
eigenvalues), but the capacity for other forms of interesting behavior might 
nevertheless remain. 

Remark 8.6 (Cyclic composition trajectories). Especially when the network 
in Corollary 18.41 is conservative, in which case stoichiometric compatibility 
classes are compact, the existence of a (positive) equilibrium that is both 
unique within its stoichiometric compatibility class and also unstable suggests 
that, for at least some some differentiably monotonic kinetics, there might 
be an attracting periodic orbit. 

There are indeed compelling reasons to believe that, for any weakly re- 
versible concordant network with discordant fully open extension, there will 

^We are grateful to Sayanti Banerjee for calling this paper to our attention. 
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exist a differentiably monotonic kinetics such the resulting kinetic system 
gives rise to a cyclic composition trajectory. In fact, certain considerations 
point to the presence of a Hopf bifurcation in connection with kinetic varia- 
tions: 

To the weakly reversible network in Corollary 18.41 one can assign a com- 
plex balanced mass action kinetics [T51 122] ■ In this case the resulting 
mass action system will have the property that the unique equilibrium in 
each nontrivial stoichiometric compatibility class has associated with it eigen- 
value^ with negative real parts. On the other hand, Theorem 18.11 ensures 
that, for the same network, there is another differentiably monotonic kinetics 
that gives rise to a nontrivial stoichiometric compatibility class containing 
a unique (positive) equilibrium associated with a positive real eigenvalue. 
Along a suitably parameterized transition from the first differentiably mono- 
tonic kinetics to the second, one therefore expects that a pair of complex- 
conjugate eigenvalues will cross the imaginary axis. (So long as the param- 
eterized kinetics remains differentiably monotonic, it follows from Theorem 
17.61 that the crossing cannot be through zero.) 

Similar considerations apply even in the case of a weakly reversible con- 
cordant network with concordant fully open extension, provided that, for the 
original network, there is a differentiably monotonic kinetics for which there 
is a positive equilibrium having eigenvalues with positive real part. However, 
for a concordant network with concordant fully open extension there is no 
guarantee, of the kind provided by Corollary 18. 4[ that such a kinetics will 
exist. 

9 Strong concordance and other generaliza- 
tions 

Motivated by a kinetic class considered in work by Banaji and Craciun, we 
introduce in this section the notion of strong concordance and examine its 
consequences. Then we consider still other variants of the concordance idea. 

^We are referring here to eigenvalues associated with eigenvectors in the stoichiometric 
subspace. 
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9.1 Strong concordance and two-way weakly mono- 
tonic kinetics 

Some — but certainly not all — of what we have done so far relied on the 
supposition of weakly monotonic kinetics, in which an increase in the rate of 
a particular reaction requires an increase in the concentration of at least one 
of its reactant species. In particular, we argued in Proposition I4.8l that. when 
the kinetics is weakly monotonic, the species-formation-rate function for a 
concordant network is invariably injective and that, as a consequence, mul- 
tiple stoichiometrically compatible equilibria, with at least one positive, are 
impossible. There was no requirement that there be a degradation reaction 
for every species. 

Although Banaji et al. [5j and Banaji and Craciun ^ do require a com- 
plete supply of degradation reactions to get injectivity results, the kinetic 
requirements they invoke are, in most ways, significantly weaker. We shall 
be more specific about the meaning of what they call nonautocatalytic (NAG) 
kinetics in Remark 19. 7[ but it suffices here to say that, with NAC kinetics, 
the rate function for a particular reaction is required to have the property 
that an increase [decrease] in the concentration of one of its reactant [prod- 
uct] species — keeping all other concentrations fixed — cannot result in a 
decrease of the reaction ratelll (Species that, for the reaction, are neither re- 
actants nor products are not permitted to have an effect on the rate.) When 
the stoichiometric matrix is "strongly sign determined" — i.e., every one 
of its square submatrices has a certain property — Banaji and Craciun are 
also able to assert injectivity and the absence of multiple steady states, but, 
again, only in the fully open setting. (In the fully open setting, the issue of 
stoichiometric compatibility becomes moot.) 

With weakly monotonic kinetics, then, the rate of a particular reaction 
can be influenced (positively) only by an increase in the concentration of a 
reactant species. The "influences" that NAC kinetics admits are broader: 
The rate of a reaction can also be influenced (positively) by a decrease in 
the concentration of a product species. In rough terms, NAC kinetics admits 
"two-way influences": It embraces product inhibition, while weakly mono- 
tonic kinetics does not. 

For the most part, it is the purpose of this section to show that, without 

^Without much loss of generality, Banaji and Craciun presume that no reaction has a 
species common to both its reactant and product complexes. 
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much difficulty, Proposition 14.81 can be generalized to subsume such two-way 
influences in the kinetics, provided that network concordance is replaced by 
strong concordance. We will assert that, when the kinetics is two-way weakly 
monotonic, the species-formation-rate function for a strongly concordant net- 
work is invariably injective and that, as a consequence, multiple stoichiomet- 
rically compatible equilibria, with at least one positive, are impossible. As 
was the case with Proposition 14.81 (and in contrast to the otherwise broad 
Banaji and Craciun results) there will be no requirement of a degradation 
reaction for every species. 

We begin by making precise what we mean by a two-way weakly mono- 
tonic kinetics. Hereafter, following Banaji and Craciun, we restrict our atten- 
tion to networks having the property that no reaction has a species common 
to both its reactant and product complexes. In practical terms, this is a very 
small price to pay for the loss of ambiguity in speaking of a reactant or 
product species for a particular reaction. 

The following deflnition is intended to describe a kinetics in which an 
increase in the rate of a particular reaction requires that there be either 
an increase in the concentration of a reactant species or a decrease in the 
concentration of a product species. And for the rate to remain unchanged, 
it is necessary that there be no change in the concentrations of any reactant 
species or else that there be appropriately opposing changes in concentrations 
of the reactant and/or product species. 

Definition 9.1. A kinetics for a reaction network {=5^,'^,^} is two-way 
weakly monotonic if, for each pair of compositions c* and c** , the following 
implications hold for each reaction y ^ y' e 1% such that supp y c supp c* and 
supp y c supp c**: 

(i) J^y^y'(c**) > J^y^yi{c*) => thcrc Is 3. specles s such that sgn (c** - c*) = 
sgn(?/-?/')s 0. 

(ii) ,Ji(fy^y'{c**) = J€y^yi{c*) => c* * = c* for all s € suppy, or else there are 
species s, s' with sgn (c** - c*) = sgn {y - y')s and sgn (c** - c*,) - 
-sgn (?/-?/')«' '^O. 

Remark 9.2. Deflnition 19 . 1 1 provides for a kinetic class that subsumes the class 
of weakly monotonic kinetics (Deflnition 14. 5p . It is intended to permit, but 
not require, the concentration of a reaction's product species to influence its 
rate. The terminology two-way weakly monotonic should not be construed to 
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suggest that there is a complete anti-symmetry with respect to the influence 
of reactant and product species: Note that — as with mass action kinetics — 
the rate of a reaction can remain unchanged even when the concentrations 
of all of the reaction's product species increases, the concentrations of all 
other species remaining fixed. The same is not true of the reaction's reactant 
species. An increase in the concentrations of all of its reactant species, the 
others remaining fixed, must result in an increase in the reaction rate. Again, 
any kinetics that is weakly monotonic (e.g., mass action kinetics) is also two- 
way weakly monotonic. 

Remark 9.3. In contrast to the definition of NAC kinetics [HIS], there is 
no requirement here that the rate functions in a two-way weakly monotonic 
kinetics be different iable or even continuous. Although Definition 19.11 is in- 
tended to capture the general spirit of NAC kinetics, there are also subtle, 
but consequential, differences that we shall discuss in Remark 19.71 

Because the class of two-way weakly monotonic kinetics is broader than 
the class of weakly monotonic kinetics, analogous theory will require focus on 
a narrower class of reaction networks — in particular the class of networks 
we call strongly concordant. In the following definition, the map L is as in 
eq. (HOD. 

Definition 9.4. A reaction network {=5^,^,^} with stoichiometric subspace 
S is strongly concordant if there do not exist a e ker L and a non-zero a e S 
having the following properties: 

(i) For each y ^ y' such that cty^y' > there exists a species s for which 
sgnas^sgn{y-y')s*0. 

(ii) For each y y' such that cty^y' < there exists a species s for which 
sgno-^ = -sgn (y-y')s 4^ 0. 

(iii) For each y y' such that cty^y' = 0, either (a) cr^ = for all s e suppy, 
or (b) there exist species s,s' for which sgnas = sgn{y - y')s + and 
sgncT^. = -sgn(y-?/')s' * 0. 

Note that a network that is strongly concordant is also concordant. Net- 
work f l25]) is concordant but not strongly concordant. 

A + B ^ P 

B + C ^ Q (28) 
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Remark 9.5 (Ascertaining strong concordance). We intend to show in a sep- 
arate article that strong concordance of a fully open network can be asserted 
when its Species-Reaction Graph satisfies the conditions stated in Remark 
12.11 More generally, strong concordance of a network, not necessarily fully 
open, can be determined, either positively or negatively, by a sign-checking 
procedure implemented in 

The following proposition generalizes Proposition 14.81 

Proposition 9.6 (Generalization of Proposition 14. 8p . A two-way weakly 
monotonic kinetic system {^,'^,^, ^} is injective whenever its underly- 
ing reaction network {^,^^,1%} is strongly concordant. In particular, if the 
underlying reaction network is strongly concordant, then the kinetic system 
cannot admit two distinct stoichiometrically compatible equilibria, at least 
one of which is positive. 

The proof of the proposition closely parallels the fairly simple proof of 
Proposition 14.81 In fact. Proposition 19.61 is merely a special case of a more 
general proposition that we state and prove in §9.2^ one that allows for very 
general species-influences in the kinetics. 

Before we do that, however, we want to emphasize once again that — 
in contrast to the work of Banaji and Cracun — Proposition 19.61 does not 
require the fully open setting. That they were compelled to require a full 
supply of degradation reactions in order to ensure injectivity of the species- 
formation-rate function is due, in part, to difficulties inherent in the NAG 
kinetic class. Although NAG kinetics is close in spirit to what we call two- 
way weakly monotonic kinetics, there are subtle distinctions between the 
two that, for NAG kinetics in the true chemical reactions, makes a result 
resembling Proposition 19.61 impossible, at least in the absence of additional 
special reactions (e.g., degradation reactions) conforming to stronger kinetic 
requirements. This we explain in the following remark. 

Remark 9.7 (Difficulties instrinsic to NAG kinetics). In the sense of [HIS], a 
(differentiable) kinetics for a reaction network {S^, ^} is non-autocatalytic 
(NAG) if for each c* € Mf^, for each reaction y y' e M, and for each species 
s € S^, the following conditions hold: 

2. y'^-y^ = 0^'-^(c*)=0. 
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For any reaction network , it is not difficult to imagine an 

NAC kinetic system J^} in which, for some open set Q e M;^, each 

rate function J^y^y'(-) is constant on Q. In this case, the species-formation- 
rate function is also constant on fl, so the kinetic system is not injective. 
Thus, there can be no assertion resembling Proposition I9.6[ in which every 
reaction network within a certain class, when taken with every NAC kinetics, 
invariably gives rise to an injective kinetic system. 

In the work of Banaji and Craciun injectivity results from the supposition 
that the "true" chemical reactions — those that conform to NAC kinetics 
— are supplemented by degradation reactions that conform to still stronger 
kinetic requirements. 

The following theorem is a variant of Theorem 17.41 It will primarily find 
use in the following way: We intend to show in a separate paper that, when a 
reaction network's Species- Reaction Graph satisfies the conditions of Remark 
12. H then the network's fully open extension is strongly concordant. Thus, 
when the network is normal, those same Species-Reaction Graph conditions 
will ensure strong concordance of the original network. 

Theorem 9.8. A normal network is strongly concordant if its fully open 
extension is strongly concordant. In particular, a weakly reversible network 
is strongly concordant if its fully open extension is strongly concordant. 

Proof of the theorem is very similar to the proof of Theorem 17.41 See 
Appendix O 

Theorem 17.61 also has a generalization. In preparation for it we need the 
following definition: 

Definition 9.9. A kinetics for a reaction network {S^ ,M} is differ- 
entiably two-way monotonic at c* e Mf^ if, for every reaction y ^ y' ^ 
^y^y'(-) is differentiable at c* and, moreover. 



(ii) ^^(c*)<0 Vs€suppy', 

(iii) ^-^^{c*)=0 Vs e =5^ \ (suppy u supp J/'). 

Proofs of the following theorem and its corollary proceed very much like 
the proofs of Theorem 17.61 and its corollary. 
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Theorem 9.10. Let be a kinetic system with stoichiomet- 

y 

ric subspace S and species formation rate function / : S. Moreover, 

suppose that the kinetics is differentiably two-way monotonic at c* e Mf". // 
the underlying network {^,^,^} is strongly concordant then the derivative 
df{c*) ■ S ^ S is nonsingular (whereupon is not one of its eigenvalues). If 
the network's fully open extension is strongly concordant then no real eigen- 
value of df{c*) is positive. 

Corollary 9.11. Let be a kinetic system with stoichiometric 

— y 

subspace S and species formation rate function / : ]R_^ S. Moreover, 
suppose that the kinetics is differentiably two-way monotonic at c* € If 
the underlying network {^,'^,^} is normal — in particular, if it is weakly 
reversible — and if the network's fully open extension is strongly concordant, 
then every real eigenvalue of df{c*) is negative. 

9.2 Still more general species influences on reaction 
rates 

In our definition of two-way weakly monotonic kinetics, we generalized weakly 
monotonic kinetics to permit (but not require) for a particular reaction an 
inhibitory effect exerted by the reaction's product species. That is, a partic- 
ular kinetics within the two-way weakly monotonic class was permitted to 
have the property that the rate of a particular reaction y ^ y' might increase 
whenever the concentration of s e supp y' decreases, the concentrations of all 
other species remaining fixed. The companion notion of strong concordance 
— an attribute of network structure alone, divorced from kinetic considera- 
tions — worked hand-in-hand with the broader kinetic class to give rise to 
the injectivity result given in Proposition 19.61 

It is not substantially more difficult to consider, in an analogous way, ki- 
netic classes with far wider ranges of permissible species influences. Indeed, 
for a particular reaction we might want to permit inducer or inhibitor effects 
on the reaction rate by species that, for the reaction, are neither reactants 
nor products. (We shall continue to require that, for each reaction, an in- 
crease in the concentration of one its reactant species, with all other species 
concentrations held constant, results in an increase of the reaction's rate.) 

We begin by indicating formally what we mean by an influence specifi- 
cation for a network. This amounts to indicating, for each reaction, which 
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species are deemed to be possible inducers and which are deemed to be possi- 
ble inhibitors (and which species are required to have no effect whatsoever). 

Definition 9.12. An influence specification J' for a reaction network "^,^1 
is an assignment to each reaction y y' of a function J^y^y' ■ y {1,0, -1} 
such that 

^y^y'(s) = 1, Vs e suppy. (29) 
If '-^y^y'{s) = 1 [resp., -1] then species s is an inducer [inhibitor] of reaction 

y^y'- 

Recall that, for the two-way monotonic kinetic class, we permitted, but 
did not require, the product species to be inhibitors for a particular kinetics 
within the class. In the following definition, our intent is to specify, for a 
broad kinetic class, which species are permitted to be inhibitors or inducers 
— without requiring them to be one or the other — for a particular kinetics 
within the class. (The exception is that, in the sense of Definition I9.12[ we 
require that, for each reaction, its reactant species are inducers.) 

Definition 9.13. A kinetics for a reaction network {o5^,^,i^}is weakly 
monotonic with respect to influence specification J" if, for every pair of 
compositions c* and c**, the following implications hold for each reaction 
y y' ^ such that supp y c supp c* and supp y c supp c**: 

(i) J^^j/'(c**) > J^y^y'ic*^ =^ there is a species s such that sgn(c** -c*) = 

(ii) ,J^y^y'(c**) = J€y^yi{c*) => clthcr (a) c** = c* for all s 6 suppy or (b) 
there are species s, s' with sgn (c**-c*) = ^y^y'{s) 4^ and sgn (c*,* - c*,) 

Remark 9.14. Note that if a kinetics Jif is weakly monotonic with respect to 
influence specification and if c** is a positive composition then implications 
(z) and {a) in Definition 19. 131 will hold for all reactions in the network, even 
when c* is not strictly positive. 

In the following definition, the map L is again as in eq. flUJ]) . 

Definition 9.15. A reaction network with stoichiometric sub- 

space 5" is concordant with respect to influence specification J' if there do not 
exist a € ker L and a non-zero o ^ S having the following properties: 
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(i) For each y y' such that dy^y' > there exists a species s for which 
sgncTs = J^y^yr{s) + 0. 

(ii) For each y ^ y' such that ay^y' < there exists a species s for which 
sgno-s = -J^y^y'(s) + 0. 

(iii) For each y ^ y' such that ay^y' = 0, either (a) cTs = for all s e supp ?y or 
(b) there are species s,s' for which sgnus = ^y^y'{s) 4^ and sgna^' = 

-^y^y,{s')4Q. 

Remark 9.16. It is expected that will be expanded to ascertain not only 
concordance and strong concordance but also the more general form of con- 
cordance indicated in Definition 19. 15[ 

Remark 9.17. The definitions of the weakly monotonic kinetic class (Defi- 
nition 14.51) and of ordinary network concordance (Definition 14.11) emerge as 
special cases of Definitions 19.131 and 19.151 applied with the influence specifi- 
cation given, for each reaction, by J'y^y' - sgny. Similarly, the definitions 
of the two-way monotonic class and of strong concordance emerge from the 
influence specification given, for each reaction, by J^y^y' = sgn {y - y'). 

The following proposition generalizes Proposition 19. 6[ In light of Remark 
I9.17[ its proof will serve also as a proof of Proposition 19. 6[ 

Proposition 9.18 (Generalization of Proposition [HS]). A kinetic system 
{^,^,^, ^} is injective whenever there exists an influence specification 
such that: 

(i) The kinetics J(f is weakly monotonic with respect to J' . 

(ii) The underlying network {^,^,^} is concordant with respect to J" . 

Proof. Let be a kinetic system and suppose that is an in- 

fluence specification that satisfies conditions (i) and (ii) in the proposition 
statement. Suppose also that, contrary to what is to be proved, there are dis- 
tinct stoichiometrically-compatible compositions c* and c**, with c** € Mf^, 
such that 

^ jey^yicn{y'-y)= E X^^yimy' -y). (30) 
Then a £ M-^, defined by 
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(31) 



is a member of kcr L. The vector a ■= c** - c* is clearly a member of the 
stoichiometric subspace, S. 

Note that, if for a particular reaction y ^ y' wc have ay^y' > 0, then weak 
monotonicity with respect to implies that, for some species s, sgno"^ - 
sgn(c** - c*) = J^y^y'{s) + 0. Similarly, when Oiy^yi < 0, there is a species s 
such that sgnus = sgn(c** - c*) = -sgn(c* - c**) = -J^y^yi{s) + 0. Finally, 
if OLy^yi = 0, then CTg = c** - c* = for all s € supp y, or else there exist 
distinct species s,s' such that sgn CTg = sgn(c** - c*) - ^y^y'{s) + and 
sgno-^/ = sgn (c^ - c^,) = -J^y^y.(s') ^ 0. 

The existence of a and cr so constructed contradicts the supposition that 
the network {5^, ^} is concordant with respect to influence specification 



It is not difficult to see that other propositions and theorems in this paper 
can also be generalized in various ways to accommodate the broader notions 
of concordance and kinetic monotonicity, both relative to a common influence 
specification. 

Remark 9.19 (Kinetic classes with mandatory inducer-inhibitor influences). 
For a reaction network , consider the class of all kinetics that are 

weakly monotonic with respect to a certain influence speciflcation, J^. More- 
over, let y be a fixed reaction, and suppose that species s is such that 
J^y^yi{s) < 0. (This choice of sign is made only for the sake of concretion.) 

Now if J(f is a fixed kinetics in the kinetic class under consideration, 
then, with respect to ^ species s might have no effect at all on the rate 
of reaction y -> y'. With another kinetics Jt' in the same class, s might 
indeed be a nontrivial inhibitor oi y y'. That J^y^y>{s) is negative merely 
indicates that the kinetic class under consideration is sufficiently wide as to 
admit some kinetics, such as J^', in which s is an inhibitor of y ^ y' . 

In this way, mass action kinetics (for which there are no inhibitors) can 
be viewed to sit within a larger class of kinetics, certain members of which 
might include inhibition. Theorems can then be stated for broad classes of 
kinetics that include mass action kinetics as a special subclass. 

We could just as well have chosen to consider kinetic classes within which 
every kinetics is required to include, for example, a nontrivial inhibition of 
species s on reaction y y'. Specification of such compulsory infiuences can 



□ 
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be accomplished in the following way: Definition 19. 121 would be expanded to 
include, for each reaction y y', a. specified species set ^y^yi c with 
supp y c ^y^yi , for which the infiuences given by J'y^y' are mandatory. 

In item (ii)(a) of Definition l9.13[ supp 2/ would then be replaced by ^y^yi^ 
and the same replacement would be made in item (iii) (a) of Definition 19.151 
With very minor adjustments in the proof, Proposition 19.181 would emerge 
as before. 
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Appendices 



A Proof of Proposition 14. 9L 



Here we prove the following proposition: 

Proposition 1331 Let {^,^,^} be a reaction network that is discordant. 
Then there exists for it a weakly monotonic kinetics Ji^ such that the resulting 
kinetic system {^,^,,^, J^} is not injective. 

Proof. Because the network {^,^,^} is discordant there exist for it a 
nonzero a in the stoichiometric subspace, S, and an a e M'* that together 
satisfy the conditions in Definition 14.11 In particular we have 

^ ay^yr{y'-y) = 0. (A.l) 

From the conditions in Definition 14.11 it is not difficult to see that, for 

— y 

each y y' e there is a vector Py^y' e such that 

supp Py^y' = supp y and Uy^yi = Py^yi ■ a. (A. 2) 
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Because sgnx = sgn (exp(x)-l) for each real x, there are positive numbers 
{Vy^y'}y^y'e.^ such that (lA.ip can be rewritten as 



Vy^y'[exp{py^y' -a) -l]{y' -y) = 0. 

For the same reason, there exists c* e Mf" such that 

as = c* [exp{as) - 1] , Vs € ^. 

We define c** e Mf by 



Then 

and 



:= c* o e". 



c - c = (T e b 



From f[0|) and f[0) it follows that 



(A.3) 

(A.4) 

(A.5) 
(A.6) 
(A.7) 

(A.8) 



Now we construct a weakly monotonic kinetics for the network in the follow- 
ing way: For each y y' e ^ we let 



'^y^y'{c):-r]y^y,{^ — ^ 

From f[01l and f[01) it follows that 

y-^y'^M y^y'^^ 

whereupon the kinetic system {^,^,^, J^} is not injective. 



(A.9) 



(A.IO) 



□ 
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B An alternative proof of Theorem 



In the proof of Theorem 16 . 31 off ered earher we required that every mass action 
system for which the underlying reaction network is weakly reversible admits 
a strictly positive equilibrium in each nontrivial stoichiometric compatibility 
class, regardless of the (positive) values that the rate constants take. That 
this is so is proved in a currently unpublished manuscript |T3]. With this in 
mind, we offer in this appendix a slightly different proof which instead draws 
on established properties of complex balanced mass action systems. (See, in 
particular, [22]. Here we draw heavily on [15j.) 
We begin with a lemma: 

Lemma B.l. Let {^,"10^,^%^} he a weakly reversible subnetwork of reac- 
tion network {^,^,^}, and let c e be a fixed but arbitrary positive 
composition. Then it is possible to find positive numbers (rate constants) 
{ky^y'}y^yi^agt tkat satisfy 

E ky^y'icny'-y)=0. (B.l) 

y~*y'eMt 

In fact, {ky^yi}y^yi^cjgt can be chosen to satisfy the stronger requirement 

E ky^y,{cy{Uy:-Uy)=Q, (B.2) 

where {uy}y^<(f is the standard basis for R*^. 

Proof. (See in particular Lecture 4 in j.l5j.) Because {^5^, ^t^i^tj is weakly 
reversible, there exist positive numbers {ny^y'}y^y'e.^t that satisfy 

Ky^y'{UJy' -Uy) = 0, (B.3) 

y^y'i^'l 

and, therefore, 

E S-.'(l/'-2/) = 0. (B.4) 
Now choose {ky^yf}y^y,^^^t to satisfy 

Ky^yr^ky^y,Cy, ^ y ^ y' ^ (B.5) 

□ 
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Remark B.2. When the rate constants have been chosen to satisfy condition 
(IB.ip . the composition c becomes an equihbrium of the mass action systems 
^'f,e^1', k}. When the rate constants satisfy the stronger condition flB.2p . 
then c is said to be an equihbrium at which complex balancing obtains, and 
{y ,k} is said to be a complex balanced mass action system. Proper- 

ties of complex balanced mass action systems are discussed in [22] and [T3] . 

With this as background we state Theorem 16.31 once again and then pro- 
vide an alternate proof. 

Theorem 16.31 For a weakly reversible concordant reaction network, no non- 
trivial stoichiometric compatibility class can have on its boundary a compo- 
sition that is reaction-transitive. 

Alternative proof of Theorem \6.!A Suppose that {^,^,=^}is a concordant 
weakly reversible reaction network and that, contrary to what is to be proved, 

y 

c* e ]R_|_ is a reaction-transitive composition on the boundary of a nontrivial 
stoichiometric compatibility class. Let 

y* := suppc*, 

'tf* := {y : supp y c supp c*} (B.6) 
'■= ^ e =^ : supp y c supp c*}. 

To be concrete, we shall begin by considering the more difficult situation 
in which ^* is not empty. By virtue of Remark l6.6[ {y* , , ^* } is a weakly 
reversible subnetwork of the original network. 

Now we choose c* to be a strictly positive composition in M.-^ that agrees 
with c* for every species in y* — that is at = c* for all s e y*. By 
virtue of Lemma IB. II and Remark IB. 21 we can choose positive rate constants 
{ky^yf}y^y'^m* to ensure that c* is a complex balanced equilibrium of the 
mass action system {^, /c*}. Because c# agrees with c* on y* and 
because every complex in ^* has support in y* we have 

E K^y'icr{y'-y)= E K^y'(c*ny'-y) = o. (b.7) 

y^y'e.^* y^y'e.'ig* 

(Note that species in \ y* are "inerts" relative to reactions in they 
are neither consumed nor produced by those reactions.) 

Because {y ,k*} is a complex balanced mass action system, it 
has the properties that accrue to such a system [T51I22]- Among these is the 
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existence of precisely one positive equilibrium in each nontrivial stoichiomet- 
ric compatibility class for the underlying network {^,'10'*,^*}. Note that 
a stoichiometric compatibility class for that network is the intersection of 

y 

with a parallel of its stoichiometric subspace, 

S* := span{y' - y e R-^ : y ^ y' € ^*}. 

On the other hand, a stoichiometric compatibility class of the parent net- 
work {^,'^,^}is the intersection of ]R_,_ with a parallel of its stoichiometric 
subspace, 

S := spanjy' -y eR^' -.y ^ y' e M}, 

Note that 5" contains S*. 

In fact, relative to , the nontrivial stoichiometric compatibility 

class (c* + S*) n containing the boundary composition c* is the union of 
stoichiometric compatibility classes for the network ,^%*}. Clearly, 

some of these contain positive compositions; those that do are nontrivial, 
and each therefore contain precisely one positive equilibrium for the mass 
action system {y ^M* ,k*}. 

Hereafter we suppose that c** is one such positive equilibrium for the 

y 

mass action system l^,"^*,^*, fc*}, residing in (c* + S) ; thus, we have 

E K^y'ic**ny'-y) = ^- (B.8) 

Having already chosen rate constants for reactions in we can now 
choose rate constants for the remaining reactions — that is, those in ^ \ 
— to ensure that c** is an equilibrium of the resulting larger mass action 
system: In fact, because the subnetwork ^ \ ^ \ is also weakly 
reversible (Remark 16 ■61) . Lemma IB.ll ensures that we can choose positive 
numbers {A;*!^y,}y^j/'e.«x,«* to satisfy 

E k;i^>{c**)y{y'-y) = Q. (B.9) 

Finally, we denote by A; e Rf the rate constant assignment for the full network 
constructed as indicated — i.e., ky^yr = k*^y, for all y y' e ^* and ky^y' = 
k** , for all y^y'€M\3g*. 
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Taken together, (IB.Sp and (IB. 90 ensure that c** is a positive equihbrium 
of the mass action system ,"1^ ,3?,k}. Moreover, c* is a (boundary) equi- 
hbrium of that same mass action system: This foUows from (IB.7P and the 
fact that, for each y ->■ y' e M \ 3?*, supp y 9^ supp c*. 

Because, relative to the concordant network {^, ^} , the distinct equi- 
hbria c* and c** are stoichiometricaUy compatible, we have a contradiction 
of Proposition I4.8[ 

We must still consider the case in which 3?* is empty. In this case, supp y ^ 
suppc* for aU y ^ so, for any kinetics assigned to {^,^,e^} , c* wih be 
an equilibrium of the resulting kinetic system. Now let c** be a positive 
composition in the nontrivial stoichiometric compatibility class containing 
c*. We can, by virtue of Lemma IB.lt assign a positive rate constant to 
each reaction of =^ in such a way as to make c** an equilibrium of the 
resulting mass action system. Then c* and c** are distinct stoichiometricaUy 
compatible equilibria of that same mass action system. Because the network 
is concordant, we again have a contradiction of Proposition 14.81 

□ 



C Proof of Theorems 17.41 and 19. 8L 



Here we provide proof of Theorem 17.41 the statement of which is repeated 
immediately below: 

Theorem 17.41 A normal network is concordant if its fully open extension 
is concordant. In particular, a weakly reversible network is concordant if its 
fully open extension is concordant. 

Proof. We suppose that the fully open extension of a normal network 'if, 
is concordant but that itself is not concordant. This will produce 

a contradiction. As usual we denote by c M-^' the stoichiometric subspace 
for the base network {^,^,^} . Note that the stoichiometric subspace for 
the fully open extension coincides with M'^. 

If the network {^,'^,^}is not concordant then, for the network, there 
is a nonzero a ^ S and an a e ker L which together satisfy the conditions in 
Definition 14.11 In particular, we have 

Y ay^y.{y'-y) = Q. (C.l) 



49 



Given such a a, a pair, our aim will be to show that the fully open extension 
{^j'tf,^} cannot be concordant, which contradicts the hypothesis. That 
is, we will want to show the existence of a nonzero a € R-^ and an a e R'^ 
satisfying 

Y, ay^y'{y'-y) = (C.2) 

such that the a, a pair also satisfies the requirements in Definition 14.11 

For the normal network {^,^,^} , let T : S* -> S* be as in Definition 17.21 
Because T is nonsingular, the map [T - e/]"^ S S exists for sufficiently 
small e. In fact, in a sufficiently small neighborhood of the function e 
[T - eiy^ is continuous, where the set of linear maps from to 5" is given 
the usual norm topology. (See [2S], §5.2.) 

For small e > 0, let x{^) ■= [T - el]~^a, which is to say that 

Tx{^) = a + ex{e). (C.3) 

(Note that x(0) = T^^a.) Thus, we can write 



E %^y'iy*x{e)W-y)+Z(^s + eXs{e)){-s)=0. (C.4) 

In anticipation of construction of a, we multiply (IC.ip by a positive num- 
ber M, to be determined later, and add the result to (]C.4p to obtain 



E {May^y,^7^y^y,{y*x{emy'-y)+Y.i^s + eXs{e)){-s)=^. (C.5) 

y^y'aSS say 

It will be useful to keep in mind that -s is the reaction vector corresponding 
to the reaction s e e;^. 
Hereafter, we choose 

a:=a + ex(e), (C.6) 

where it is understood that e > is chosen sufficiently small as to satisfy the 
following requirement: 

sgn (cr + ex(e))s = sgnds when sgnas 0. (C.7) 

Because a i^Q, flC.7p ensures that a 0. 
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We turn now to a choice of a. By way of preparation, note that we 
have not precluded the possibihty that ^ might contain species-degradation 
reactions of the form s ^ 0, where s is a species. In fact, let ^ c ^ be the 
set of all species for which there are such reactions; that is, = {s e : 
s ^ € We denote by the complement of ^ in J^. By we mean 
the set of all reactions in M that are not species-degradation reactions; that 
is, is the set of all members of ^ not of the form s ^ 0, s e ^ . 

With this as background, we choose a e R-^ as follows: 

ay^y' := May^y' +'r]y^yr{y * xi(^)) ^oi dXl ^ y' ^ , (C.8) 
as^Q ■■= Mas^o + Vs^ois * x(e)) + (a.^ + eXs(e)) for all s € 
«s^o ■= (o"s + eXs(e)) for all s 6 , 

and we suppose that M > is sufficiently large as to satisfy the following 
requirement: For all y ^ y' e M 

sgnay^y' = sgn ay^y' if sgn ay^y' 4^ 0. (C.9) 

Note that, by virtue of (IC.Sp . a so chosen is a solution of (]C.2p . 

It remains to be shown that this a, a pair, constructed as indicated, sat- 
isfies the conditions in Definition I4.1[ There are several cases to consider, 
corresponding to various categories of reactions in 

1. For reactions of the form s -> 0,s e we have sgncis^o = sgna^ so, for 
them, the conditions in Definition 14.11 are clearly satisfied. 

2. For reactions of the form s ^ 0, s € there are two possibilities: 

(a) If as^Q then we have sgna^^o = sgna^^o - sgncXs = sgna^ so 
Condition (i) in Definition 14.11 is satisfied. 

(b) If as^o = then 0-^=0 and sgna^^o = sgn (ijs^ois * x(e)) + exs{,e)) = 
sgn (e) = sgn (ct^ + e^s (e) ) = sgn a, . 

Therefore, either Condition (i) or Condition (ii) in Definition 14.11 is satis- 
fied. 

3. For y -y y' ^ such that ay^y' there must exist s e supp y such that 
sgn ay^y' = sgn ay-^y' = sgncr^ = sgno-g, so that Condition (i) in Definition 
14. II is satisfied. 

4. U y y' e is such that cty^y' = and supp y n supp a = then 

sgnoy^y. = sgn(y*x(e)), (C.IO) 
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and, for all s e supp y, 



sgna^ = sgn(as + exs(e)) =sgnXs(e)- (CH) 
There are two cases to consider: 



a) li y * then, from ( IC.lOp . we have the existence of s' e suppy 



such that sgna^^y. = sgny * x(e) = sgnxs'(e) = sgn{as' + exs'i^)) = sgncr^ 
Thus, Condition (i) in Definition 14.11 is satisfied. 

(b) If y * = then either suppi/ n suppx(e) = 0, in which case 
supp y n supp a = 0, or there exist s, s' € supp y such that sgn Xs(,^) = -sgn Xs'(e), 
both nonzero, in which case we have sgn ag = -sgn dg', both nonzero. In either 
case, then. Condition (ii) in Definition 14.11 is satisfied. 

5. Again consider y y' e 3?' such that ay^yi = 0, but now suppose and 
there exist s,s' € supp?/ such that sgncTs = -sgncs', both non-zero. Then 
from flC.7p we have sgna^ = -sgna^/, both non-zero. Regardless of the sign 



of Ciy^y', therefore, the conditions in Definition 14 . 1 1 are satisfied. 

In summary, then, the a, a pair, constructed as indicated, ensures that 
l^,^,,^}, the fully open extension of the (presumed discordant) normal 
network {S^ is not concordant. This, however, contradicts what 
has been supposed about the fully open extension. Therefore, {^,^,=^}is 
concordant. □ 

We repeat the statement of Theorem 19.81 below: 

Theorem 19.81 A normal network is strongly concordant if its fully open 
extension is strongly concordant. In particular, a weakly reversible network 
is strongly concordant if its fully open extension is strongly concordant. 

Proof of the theorem is very similar to the proof of Theorem 17. 4[ The 
constructions of a and a are identical. The only substantive changes are in 
items 3 and 5: 

3. For y ^ y' e such that ay^y' 4^ there must exist a species s such that 
4^ {sgn Oy^yr) {sgn (y - y')s) = {sgn ay^y,) {sgn {y - y')s) = sgna, = sgna^, so 
that Conditions (i) and (ii) in Definition 19.41 are satisfied. 

5. Consider y ^ y' e 3?' such that cty^y' = 0, but now suppose and there exist 
species s,s' such that 

sgno-^ =sgn{y-y')s 4 
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sgaa.r = -sgia{y -y')s' 4^ 0. 

From (1C.7P we also have sgnas = sgncTs and sgnas' = sgncis'. Regardless of 
the sign of dy^y', then, the conditions in Definition 19.41 are satisfied. 



D Proofs of Theorems 7.6 and 8.1 



We repeat below the statement of Theorem 17. 6t 

Theorem 17.61 Let be a kinetic system with stoichiometric 

y 

subspace S and species formation rate function / : ]R_,_ -> 5*. Moreover, sup- 
pose that the kinetics is differentiably monotonic at c* £ Mf^. If the underlying 
network {y^'^^M^is concordant then the derivative df(c*) : S ^ S is non- 
singular (whereupon is not one of its eigenvalues). If the network's fully 
open extension is concordant then no real eigenvalue of df(c*) is positive. 

^ 

Proof. For the sake of brevity, for each y y' e ^ we denote by Py^y' e M,^ 
the gradient of J^^j,'(-) at c*. That is, for each s € 

K^, (D.i) 

For each y ^ y' ^ M, note that supp Py^y' = suppy. In terms of the notation 
we have now introduced, df{c*) takes the following form: For each a ^ S, 

df{c*)a= ^ {py^y,-a){y'-y). (D.2) 

Suppose that {=5^,^,^} is concordant. Suppose also, contrary to what 
is to be proved, that df(c*) is singular. Then there is a nonzero a ^ S such 
that 

E {Py^ycr){y'-y) = 0. (D.3) 

If we take a e M-^ to be defined by ay^y' := Py^y' ■ cr, then it is not difficult 
to see that, taken together, a and a satisfy the conditions of Definition 14.11 
Thus, {^,'^,^}is discordant, and we have a contradiction. 

Now suppose that the fully open extension of is concordant 

and that, contrary to what is to be proved, df(c*) has a positive real eigen- 
value. Then there is a nonzero a e S and a positive real number A such that 
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df(c*)a = Act, which can be written in the following way: 

E iPy^y^)iy'-y)+i:^^s{-s)=0. (D.4) 

Note that -s is the reaction vector for reaction s in the fully open 
extension {^,^,^} of the network {^,^,^} . (Recall the notation given 
in ([22D.) 

It is possible that the reaction set for the original network, contains 
species-degradation reactions of the form s 0, where s is a species. Let 
^ c ^ be the set of all species for which there are such reactions; that is, 
^^{sey-. s-^Oe We denote by the complement of ^ in 
By M' we mean the set of all reactions in M that are not species-degradation 
reactions; that is, is the set of all members of ^ that are not of the form 
s ^ 0, s € Then flD.4p takes the form 



E (Py^y' ■ ^)(y' - y) + E (P^-o + A)a,(-s) + E Aa,(-s) = 0. (D.5) 

Now we choose a e M-* to be defined by cty^y' := Py^y' • 5" for y ^ y' 
o^s^o '■= {Ps^o + A)5"s for s € and Os^o •= ^5"^ for s e . Thus, a satisfies 

E ay^y'(y'-y) = 0. 

It is easy to see that, taken together, b and a satisfy the conditions of 
Definition 14.11 Thus, the fully open extension {^,^,e^} is discordant, and 
we again have a contradiction. It follows, then, that no real eigenvalue of 
df{c*) is positive. □ 

We repeat below the statement of Theorem I8.lt 

Theorem 18.11 Consider a reaction network with positively dependent reac- 
tion vectors. If the network is discordant, then there exists for it a differ- 
entiably monotonic kinetics such that the resulting kinetic system admits a 
positive degenerate equilibrium. If the network's fully open extension is dis- 
cordant then there exists for the original network a differentiahly monotonic 
kinetics such that the resulting kinetic system admits an unstable positive 
equilibrium — in fact, a positive equilibrium associated with a positive real 
eigenvalue. 



54 



Proof. Throughout this proof we consider a reaction network ,M} 
with stoichiometric subspace S, and we suppose that the reaction vectors are 
positively dependent. Thus, we have positive numbers {Ky^y'}y^y'eM that 
satisfy 

^ Ky^y,{y'-y) = 0. (D.6) 
We divide the proof into two parts: 

(i) First we suppose that the network is itself discordant. 

Then there exists a pair, a € M-* and nonzero a e S, that together satisfy 

conditions (i) and (ii) in Definition 14.11 From those conditions it is not 

— y 

difficult to see that, for each reaction y ^ y' e ^ there is a vector Py^y' e 
such that 

supppj^^j^' = supp y and ay^y' = Py^y' ■ a. (D.7) 
Thus we can write 

y 

Let c* be any member of Mf^, and for each y ^ y' e ^ let Qy^y' e be 
defined by 

Furthermore let k e be defined by the requirement that 

ky^y,{c*yy-y' = Ky^y, V?/ ^ € ^. (D.IO) 

We construct a differentiably monotonic kinetics for the network as follows: 
For each y ^ y' e ^ 

J^y^y' {c) := ky^y, diy-y' , Vcelf. (D-H) 
With this kinetics the species-formation-rate function takes the form 

/(c) := E ky^y,c'^y-y'{y'-y). (D.12) 

From f lDlej) and (Imoj) it follows that, for the kinetic system {=5^,"^,^,^}, 
c* is a positive equilibrium. Taken with (ID.lOp and f lD.9|) . some computation 
gives 

df{c*)a= {Py^y<y)W-y)- (D-13) 
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Because a is not zero, it follows from (ID.Sp that df{c*) is singular, whereupon 
c* is a degenerate positive equilibrium. 

(ii) Next we suppose instead that the fully open extension of {y, ^} is 
discordant. Then there are numbers {a^^y'ly^^'e,* and {as^o}se.y satisfying 

Y, ay^yiy' - y) + ^s^o{-s) = (D.14) 

and also a nonzero a such that conditions (i) and (ii) of Definition 14.11 are 
satisfiecl§. Again we can choose for each reaction y ^ y' e M a, vector Py^y' e 

y 

]R_^ such that 

supp Py^y' = supp y and ay^y' = Py^y' ■ a. (D.15) 
Moreover, we can choose positive numbers {ms}s^y such that 

as^o = msas. (D.16) 

Thus, (1D.14P can be rewritten 

E (Py^y'-(^)(y' -y)= Jl^^^^(^)- 

Now let 

a:=Em,a,(s). (D.18) 

From flD.17p if follows that a is a member of the stoichiometric subspace 
for the original network . Moreover, because o" 0, it is easy to 

see that cr • a > 0, whereupon a i^Q. 

_ y 

If, for each reaction y y' e we let Py^y' e be defined by 

fy^y'--=^, Vse^, (D.19) 

then (ID. 171) can be rewritten as 

y^y'e.^ 



^We have not precluded the possibihty that a degradation reaction, say s 0, might 
be a member of the set of reactions in the original network. In this case the term 
as^o{-s) would appear twice in (|D.14p . once in the first sum and once in the second. 
Nevertheless, it is easy to confirm that the sentence remains true as written. 
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Hereafter, this part of the proof will very much resemble part (i): Let c* 
be any member of ]R:f , and for each y ^ y' lei Qy^y' e be defined by 

As before, we let k e M^f be defined by the requirement that 

ky^y,{c*Yy-y' = Ky^y,, ^ y ^ y' ^ (D.22) 

and also as before we construct a differentiably monotonic kinetics for the 
network in the following way: For each y^y'^^ 

jey^y> {c) -.^ ky^y, d^y^y' , M c . (D.23) 
Again, the species-formation-rate function takes the form 

/(c) ^ ky^y,diy-y'{y'-y). (D.24) 

Note that c* is an equihbrium of the kinetic systems {^,^,^, ^}, and 
some computation gives 

df{c*)a = a. (D.25) 

Thus, a is an eigenvector of df{c*) corresponding to the positive real eigen- 
value 1. Thus, the equilibrium c* is unstable. □ 
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